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ABSTRACT 

We present a catalog of high-energy gamma-ray sources detected by the Large 
Area Telescope (LAT), the primary science instrument on the Fermi Gamma-ray 
Space Telescope (Fermi), during the first f f months of the science phase of the 
mission, which began on 2008 August 4. The First Fermi-L AT catalog (IFGL) 
contains 1451 sources detected and characterized in the 100 MeV to 100 GeV 
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range. Source detection was based on the average flux over the 11-month period, 
and the threshold likelihood Test Statistic is 25, corresponding to a significance 
of just over 4a. The IFGL catalog includes source location regions, defined in 
terms of elliptical fits to the 95% confidence regions and power-law spectral fits 
as well as fiux measurements in 5 energy bands for each source. In addition, 
monthly light curves are provided. Using a protocol defined before launch we 
have tested for several populations of gamma-ray sources among the sources in the 
catalog. For individual LAT-detected sources we provide firm identifications or 
plausible associations with sources in other astronomical catalogs. Identifications 
are based on correlated variability with counterparts at other wavelengths, or on 
spin or orbital periodicity. For the catalogs and association criteria that we have 
selected, 630 of the sources are unassociated. Care was taken to characterize the 
sensitivity of the results to the model of interstellar diffuse gamma-ray emission 
used to model the bright foreground, with the result that 161 sources at low 
Galactic latitudes and toward bright local interstellar clouds are flagged as having 
properties that are strongly dependent on the model or as potentially being due 
to incorrectly modeled structure in the Galactic diffuse emission. 

Subject headings: Gamma rays: observations — surveys — catalogs; Fermi 
Gamma-ray Space Telescope; PACS: 95.85.Pw, 98.70.Rz 

1. Introduction 

The Fermi Gamma-Ray Space Telescope has been routinely surveying the sky with the 
Large Area Telescope (LAT) since the science phase of the mission began in 2008 August. 
The combination of deep and fairly uniform exposure, good per-photon angular resolution, 
and stable response of the LAT have made for the most sensitive, best-resolved survey of 
the sky to date in the 100 MeV to 100 GeV energy range. 

Observations at these high energies reveal non-thermal sources and a wide range of 
processes by which Nature accelerates particles. The utility of a uniformly-analyzed catalog 
such as this is both for identifying special sources of interest for further study and for 
characterizing populations of 7-ray emitters. The LAT survey data analyzed here allow much 
more detailed characterizations of variability and spectral shapes than has been possible 
before. 



Here we expand on the Bright Source List ( lAbdo et al.l l2009nl . BSL), which was an 
early release of 205 high-signiflcance (likelihood Test Statistic TS >100; see § 14. 3p sources 
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detected with the first 3 months of science data. The expansion is in terms of time interval 
considered (11 months vs. 3 months), energy range (100 MeV - 100 GeV vs. 200 MeV - 
100 GeV), significance threshold {TS > 25 vs. TS > 100), and detail provided for each 
source. Regarding the latter, we provide elliptical fits to the confidence regions for source 
location (vs. radii of circular approximations), fiuxes in 5 bands (vs. 2 for the BSL) for the 
range 100 MeV - 100 GeV, and monthly light curves for the integral fiux over that range. 



We also provide associations with previous 7-ray catalogs, for EGRET (IHartman et al. 



19991 : ICasandiian fc Grenierll2008r i and AGILE (IPittori et al.ll2009f ). and with likely counter- 
part sources from known or suspected source classes. The number of sources for which no 
plausible associations are found is 630, at the spe cified confidence le vel for source associa- 
tion (80%). The First LAT AGN Catalog (ILAC, lAbdo et al.lboiOlh is based on the IFGL 
sources, and applies the same association methods, but provides associations for AGNs down 
to the 50% confidence level. 

As with the BSL, the First Fermi-LAT catalog of 7-ray sources (IFGL, for first Fermi 
Gamma-ray LAT) is not fiux limited and hence not uniform. As described in § HI the 
sensitivity limit depends on the region of the sky and on the hardness of the spectrum. Only 
sources with TS > 25 (corresponding to just over 4 a statistical significance) are included, 
as described below. 



2. Gamma-ray Detection with the Large Area Telescope 



The LAT is a pair-production telescope ( lAtwood et al.l 120091 ). The tracking section has 
36 layers of silicon strip detectors to record the tracks of charged particles, interleaved with 
16 layers of tungsten foil (12 thin layers, 0.03 radiation length, at the top or Front of the 
instrument, followed by 4 thick layers, 0.18 radiation length, in the Back section) to promote 
7-ray pair conversion. Beneath the tracker is a calorimeter comprised of an 8-layer array 
of Csl crystals (1.08 radiation length per layer) to determine the 7-ray energy. The tracker 
is surrounded by segmented charged-particle anticoincidence detectors (plastic scintillators 
with photomultiplier tubes) to reject cosmic-ray background events. The LAT's improved 
sensitivity compared to EGRET stems from a large peak effective area (~8000 cm^, or 
~6 times greater than EGRET's), large field of view (~2.4 sr, or nearly 5 times greater 
than egret's), good background rejection, superior angular resolution (68% containment 
angle ~0.6° at 1 GeV for the Front secti on and about a factor of 2 larger for the Back 
section, vs. ~1.7° at 1 GeV for EGRET; [Thompson et al.lll993l ). and improved observing 
efficiency (keeping the sky in the field of view with scanning observations, vs. inertial 
pointing for EGRET). Pre-launch predictions of the instrument performance are described 
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in lAtwood et all f l2009h . 



The data analyzed for the IFGL catalog were obtained during 2008 August 4 - 2009 
July 4 (LAT runs 239557414 through 268411953, where the numbers refer to the Mission 
Elapsed Time (MET) in seconds since 00:00 UTC on 1 January 2001, at the start of the data 
acquisition runs). During most of this time Fermi was operated in sky-scanning survey mode 
(viewing direction rocking 35° north and south of the zenith on alternate orbits). During 
May 7-20 the rocking angle was increased to 39° for operational reasons. In addition, a few 
hours of special calibration observations during which the rocking angle was much larger 
than nominal for survey mode or the configuration of the LAT was different from normal 
for science operations were obtained during the period analyzed. Time intervals when the 
rocking angle was larger than 43° have been excluded from the analysis, because the bright 
limb of the Earth enters the field of view (see below). 

In addition, two short time intervals associated with 7-ray bursts (GRB) that were 
detected in the LAT hav e been excluded. T hese intervals correspond to GRB 080916C (MET 



243216749-243217 979, lAbdo et all 120091 1 and GRB 090510 (MET 263607771-263625987, 



Abdo et al.ll2009a[ ). 



Observations were nearly continuous during the survey interval, although a few data 
gaps are present due to operational issues, special calibration runs, or in rare cases, data loss 
in transmission. Table [T] lists all data gaps longer than 1 h. The longest gap by far is 3.9 d 
starting early on March 16; together the gaps longer than 1 h amount to ~7.9 d or 2.4% of 
the interval analyzed for the IFGL Catalog. 

The total live time included is 245.6 days (21.22 Ms). This corresponds to an absolute 
efficiency of 73.5%. Most of the inefficiency is due to time lost during passages through the 
South Atlantic Anomaly (~13%) and to readout dead time (9.2%). 

The standard onboard fil tering, event reconstruction, and classification were applied to 
the data ( Atwood et al. 2009), and for this analysis the 'Diffuse' event clas^ is used. This 



is the class with the least residual contamination from charged-particle background events, 
released to the public. The tradeoff for using this event class relative to the 'looser' Source 
class is primarily reduced effective area, especially below 500 MeV. 

The instrument response functions (IRFs) - effective area, energy redistribution, and 
point-spread function (PSF) - used in the likelihood analyses described below were derived 
from GEANT4-based Monte Carlo simulations of the LAT using the event-selection cri- 
teria corresponding to the Diffuse event class. The Monte Carlo simulations themselves 



See http://fermi.gsfc.nasa.gov/ssc/data/analysis/docuinentation/Cicerone/CiceroneJData/LAT_DP.html 
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Table 1. Gaps Longer Than One Hour in Data 



Start of Gap (UTC) Duration (h) 



2008-09-30 


14:27 


1.16 


2008-10-11 


03:14 


1.59 


2008-10-11 


11:04 


3.82 


2008-10-14 


12:23 


3.83 


2008-10-14 


17:11 


3.49 


2008-10-14 


20:22 


1.59 


2008-10-15 


17:03 


3.47 


2008-10-16 


15:18 


1.83 


2008-10-22 


19:20 


1.59 


2008-10-22 


23:43 


2.06 


2008-10-23 


11:16 


1.91 


2008-10-30 


16:43 


1.59 


2008-12-11 


17:41 


6.37 


2009-01-01 


00:35 


1.72 


2009-01-06 


20:43 


6.98 


2009-01-13 


13:26 


2.10 


2009-01-17 


12:58 


2.05 


2009-01-28 


19:28 


4.78 


2009-02-01 


15:46 


1.59 


2009-02-15 


10:15 


1.05 


2009-03-16 


00:27 


116.78 


2009-05-02 


19:04 


8.94 


2009-05-07 


15:21 


5.46 


2009-06-26 


12:59 


3.19 
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were calibrated prio r to launch using accelerator tests of flight-spare 'towers' of the LAT 
( lAtwood et al.l l2009l ) and have since been up dated based on ob servation of pile-up effects 
on the reconstruction efficiency in flight data (iRando et al.ll2009[ ). The effect introduces an 
inefficiency that is proportional to the trigger rate and dependent on energy. The likelihood 
analysis for characterizing the sources uses the P6_V3 IRFs (see § 14.31) . which have the ef- 
fective areas corrected for the inefficiency corresponding to the overall average trigger rate 
seen by the LAT. The use of the P6_V3 IRFs allows the energy range of the analysis for 
the catalog to be extended down to 100 MeV (vs. 200 MeV for the BSL analysis, which 
used P6_V1). Below 100 MeV the effective area is relatively small and strongly dependent 
on energy. These considerations, together with the increasing breadth of the PSF at low 
energies (scaling approximately as 0.8°(i?/lGeV)~°'^), motivated the selection of 100 MeV 
as the lower limit for this analysis. 



The alignment of the Fermi observatory viewing dire ction with the z-ax is of the LAT 
was found to be stable during survey- mode observations (lAbdo et al.ll2009rl ). Analyses of 
flight data suggest that the PSF is somewhat broader than the calculated Diffuse class PSF 
at energies greater than ~10 GeV; the primary effect on the current analysis is to decrease 
the localization capability somewhat. As discussed below, this is taken into account in the 
catalog by increasing the derived sizes of source location regions by 10%. 

For the analysis, a cut on zenith angle (angle between the boresight of the LAT and the 
local zenith) was applied to the Diffuse class events to limit the contamination from albedo 
7-rays from interactions of cosmic rays with the upper atmosphere of the Earth. These 
interactions make the limb of the Earth (zenith a ngle ~113° at the 565 km, nearly-circular 
orbit of Fermi) an intensely-bright 7-ray source (IThompson et al.lll98ll ). The limb is very 
far off axis in survey-mode observations, at least 70° for the data set considered here because 
of the rocking angle requirement described above. Removing events at zenith angles greater 
than 105° affects the exposure calculation negligibly but reduces the overall background rate. 
After these cuts, the data set contains 1.1 x 10"^ Diffuse-class events with energies >100 MeV. 

The intensity map of Figure [1] summarizes the data set used for this analysis and shows 
the dramatic increase of the brightness of the 7-ray sky at low Galactic latitudes. The 
corresponding exposure is relatively flat and featureless as was the case for the shorter time 
interval analyzed for the BSL. The degree of exposure nonuniformity is relatively small 
(about 30% difference between minimum and maximum), with the deflcit around the south 
celestial p ole due to loss of ex posure during passages of Fermi through the South Atlantic 
Anomaly JAtwood et al.1120091 ). 



3. Diffuse emission model 



An essential input to the analyses for detecting and characterizing 7-ray sources in 
the LAT data is a model of the diffuse 7-ray intensity of the sky. Interactions between 
cosmic rays and interstellar gas and photons make the Milky Way a bright, structured 
celestial foreground. Unresolved emission from extragalactic sources contributes an isotropic 
component as well. In addition, residual charged-particle background, i.e., cosmic rays that 
trigger the LAT and are misclassified as 7-rays, provides another approximately isotropic 
background. For the analyses described in this paper we used models for the Galactic diffuse 
emission (gll_iem_v02 . f it) and isotropic backgrounds that were developed by the LAT 
team and made publicly available as models recommended for high-level analyses. The 
models, along with descriptions of their derivation, are available from the Fermi Science 
Support Centei0. 

Briefly, the model for the Galactic diffuse emission was developed using spectral line 
surveys of H I and CO (as a tracer of H2) to derive the distribution of interstellar gas in 
Galactocentric rings. Infrared tracers of dust column density were used to correct column 
densities as needed, e.g., in directions where the optical depth of H I was either over or 
under-estimated. The model of the diffuse 7-ray emission was then constructed by fitting 
the 7-ray emissivities of the rings in several energy bands to the LAT observations. The 
fitting also required a model of the inverse Compton emission that was calculated using 



GALPROP (IStrong et al.ll2004j : IStrona 120071 ) and a model for the isotropic diffuse emission. 



The isotropic component was derived as the residual of a fit of the Galactic diffuse 
emission model to the LAT data at Galactic latitudes above |6| = 30° and so by construction 
includes the contribution of residual (misclassified) cosmic rays for the event analysis class 
used (Pass 6 Diffuse; see § [2]). Treating the residual charged particles as effectively an 
isotropic component of the 7-ray sky brightness rests on the assumption that the acceptance 
for residual cosmic rays is the same as for 7-rays. This approximation has been found to be 
acceptable; the numbers of residual cosmic-ray background events scale as the overall livetime 
and any acceptance differences from 7-rays would not introduce small-scale structure in the 
models for likelihood analysis. 



^http://fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels. html 



- 12 - 



4. Construction of the Catalog 



The procedure used to build the IFGL catalog follows the same steps described in 
Abdo et al.l (l2009nl ) for the BSL, with a number of improvements. We review those steps in 
this section, highlighting what was done differently for IFGL. 

Three steps were applied in sequence: detection, localization, significance estimation. 
In this scheme the threshold for inclusion in IFGL is defined at the last step, but the 
completeness is controlled by the first one. After the list was defined we determined the 
source characteristics (flux in 5 energy bands, time variability). The IFGL catalog includes 
much more information for each source than the BSL. In what follows, flux F means photon 
flux and spectral index F is for photons (i.e., F oc E"^). 

In constructing the catalog the source detection step was applied only to the data from 
the full 11-month period as a whole. That is to say, we did not search for potentially flaring 
sources which might only be detectable o n shorter timescales . Independently of this work, 
the LAT Automated Science Processing (lAtwood et al.ll2009l ) and Flare Advocate activity 
provide a framework through which such flaring sources are detected in a timely manner and 
reported as Astronomer's Telegrams (ATels). However, since all bright flaring sources that 
were reported as ATels were also bright enough to be detected over 11 months based on their 
average fluxes, they are included in the IFGL catalog anyway. No GRB is detected over the 
full interva; the time ranges of the two brightest GRBs were excluded from the analysis (see 
§[2l)- 



The pulsars (lAbdo et al.ll2010m[ ) and X-ray binaries (I Abdo et al.ll2009kl JdJql) which are 
identifled via their rotation or orbital period, were detected and localized as ordinary sources. 
But they were entered explicitly at their true positions in the main maximum likelihood 
analysis (§ 14. 3p . in order not to bias their characteristics and those of their surroundings if 
the Galactic diffuse model is impe rfect [3l). For the LAT-detected pulsars, we used the 
radio or 7-ray timing localization (lAbdo et al.l 1201 0ml ) which is always more precise than 
that based on the spatial distribution of the events. We have checked that the positions 
of the brightest pulsars found by the localization algorithm (§ 14. 2p were consistent with 
their true positions at the 95% level (using only the statistical error, without any systematic 
correction). 



4.1. Detection 



The detection step used the same ideas that were detailed in I Abdo et al.l (l2009n( ). It 
was based on the same three energy bands, combining Front and Back events to preserve 
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spatial resolution. The detection does not use events below 200 MeV, which have poor 
angular resolution. It uses events up to 100 GeV. The full band (6.7 x 10^ counts) starts at 
200 MeV for Front and 400 MeV for Back events. The medium band (12.0 x 10^ counts) 
starts at 1 GeV for Front and 2 GeV for Back events. The hard band (10.7 x 10^ counts) 
starts at 5 GeV for Front and 10 GeV for Back events. 

We used the same partitioning of the sky into 24 planar projections as in the BSL, 
and the s ame two wavelet-based detection metho ds: mr_filter ( IStarck fc Pierrd 119981 ) and 
PGWave (jPamiani et al.l 119971 : ICiprini et al.l 120071 1. The methods looked for sources on top 
of the diffuse emission model described in § [31 For mr-fiUer the threshold w as set in each 
image using the False Discovery Rate procedure (IBenjamini fc Hochberg||l995l ) at 5% of false 
detections. For PGWave we used a flat threshold at 4 a. For comparison with the BSL, 
the number of 'seed' sources from mr_filter was 857 in the full band, 932 in the medium 
band and 331 in the hard band. Contrary to the BSL procedure, we combined the results 
of those two methods (eliminating duplicates) rather than choosing a baseline method and 
using the other for comparison. The rationale was to limit the number of missed sources to 
a minimum, since the later steps do not introduce any additional sources. Duplicates were 
defined after the first localization {pointfit in § 14. 2[ run separately on each list of seeds). If 
two resulting positions were consistent within the quadratic sum of 95% error radii only one 
source was kept (that with highest significance estimate). Where pointfit did not converge, 
the 95% error radius was set to 0.3°, typical for faint sources (§ 14. 2p . 

To that same end we also introduced for IFGL two other detection methods: 



pointfind, a tool that searches for candidate point sources by maximi zing the like- 



lihoo d function for trial point sources at each direction in a HEALPix ( IGorski et al. 
20051 ) order 9 (pixel size ~0.1 deg^) tessellation of the sky. The algorithm for evalu- 



ating the likelihood is optimized for speed by using energy- dependent binning of the 
photon data, choosing 4 energy bands per decade starting at 700 MeV, and a HEALPix 
order commensurate with the PSF width in each band. A first pass examines the sig- 
nificance of a trial point source at the center of each pixel, on the assumption that 
the diffuse background is adequately described by the model for Galactic diffuse emis- 
sion and ignoring any nearby point sources. The likelihood is optimized with respect 
to the signal fraction (i.e., the source and diffuse intensities are not fit separately) 
in each energy band, with the total likelihood being the product over all the bands. 
This makes the result independent of the spectrum of the point source or of the diffuse 
background. While the search is quite efficient, it produces many false signals, so a sec- 
ond pass is used to optimize a more detailed likelihood function which includes nearby 
detected sources and fits the test source flux and diffuse background normalization 
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independently. The result of the second pass is a map of Test Statistic from which the 
coordinates of candidate point sources can be derived. 



the minimum spanning tree (jCampana et al.l |2008[ ) looks for clusters of high-energy 
events (> 4 GeV outside the Galactic plane and > 10 GeV at |6| < 15°). It is re- 
stricted to high energies because it does not account for structured background, but 
can efficiently detect very hard sources. 



We combined the 'seed' positions from those two methods with those from the wavelet-based 
methods, using the same procedure for removing duplicates as above. 

Finally, we intr oduced external seeds from the BZCAT (IMassaro et al. I boooi l and WMAP 
( IWright et al.ll2009l ) catalogs. The BZCAT catalog is not homogenous but inclu des the great 



majo rity of known, well-characterized blazars. It is a superset of the CGRaBS (IHealey et al. 



20081 ) catalog and has broader sky coverage. The WMAP catalog includes mainly bright 
FSRQ blazars, and was used primarily to try to recover soft-spectrum sources that might 
have been missed by the source-detection algorithms. 

In order to not bias the IFGL catalog toward those external sources, we used them as 
seeds only when there was no seed from the detection methods within its 95% error radius. 
Of the 335 BZCAT seeds introduced, 24 survived as LAT 7-ray sources in this catalog. Of 
the 7 WMAP seeds, 3 remain in the catalog. 

The variety of seeds that we used means that the catalog is not homogeneous. Because 
of the strong underlying diffuse emission, achieving a truly homogeneous catalog was not 
possible in any case. Our aim was to provide enough seeds to allow the main maximum 
likelihood analysis (§ 14. 3p to be the defining step of the catalog construction. The total 
number of seeds was 2433. 



4.2. Localization 

The localization of faint or soft sources is more sensitive to the diffuse emission and to 
nearby sources than for brighter sources, so we proceeded in three steps instead of just one 
for the bright sources considered in the BSL: 

1. The first step consisted of localizing the sources before the main maximum likelihood 
analysis (§ 14. 3p as we did for the BSL (using pointfit), treating each source indepen- 
dently but in descending order of significance and incorporating the bright sources into 
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the background for the fainter ones. This is fast and provides a good enough starting 
point for step 2. 

2. The second step consisted of improving the locahzation within the main maximum 
hkehhood analysis (§ I4.3P using the gtfindsrc utihty in the Science Tool^. Again 
sources are considered in descending order of significance. When locahzing one source, 
the others are fixed in position, but the fluxes and spectral indicies of sources within 
2° are left free to accommodate the loss of low energy photons in the model if the 
source that is being localized moves away. At the end of that step we have a good 
representation of the location, flux and spectral shape of the sources over the entire 
sky, but a single error radius to describe the error box. 

3. The third step is new and described in more detail below. It uses a similar framework 
as the first step, but incorporates the results of the main maximum likelihood analysis 
for all sources other than the one being considered, so it has a good representation of 
the source's surroundings. It is faster than gtfindsrc and gttsmap and returns a full 
Test Statistic map around each source and an elliptical representation as well as an 
indicator of the quality of the elliptical fit. 



The first and third steps used a hkelihood analysis tool {pointfit) that provides speed at 
little sacrifice of precision by maximizing a specially-constructed binned likelihood function. 
Photons are assigned to twelve energy bands (four per decade from 100 MeV to 100 GeV) 
and HEALpix-based spatial bins for which the size is selected to be small compared with the 
scale set by the PSF. Since the PSF for Front-converting photons is significantly smaller 
than that for Back conversions, there are separate spatial bins for Front and Back. Note 
that the width of the PSF at a given energy is only a weak function of incidence angle. For 
pointfit the likelihood function is evaluated using the PSF averaged over the full field of view 
for each energy band. For each band, we define the likelihood as a function of the position 
and flux of the assumed point source, and adopt as the background the sum of Galactic 
diffuse, isotropic diffuse (see §|3]) and any nearby (i.e., within 5°), other point sources in the 
catalog. The flux for each band is then evaluated by maximizing the likelihood of the data 
given the model using the coordinates defined by gtfindsrc. The overall likelihood function, 
as a function of the source position, is then the product of the band likelihoods. We define 
a function of the position p, as 2(log(Lmax) — log(-^(p)), wh ere L is the likelihood function 
described above. This function, according to Wilks' theorem (jWilkslll938l ). is the probability 
distribution for the coordinates of the point source consistent with the observed data. Note 



•^Available from the Fermi Science Support Center, 



http://fermi.gsfc.nasa.gov/ssc 
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that the width of this distribution is a measure of the uncertainty, and that it scales directly 
with the width of the PSF. 

We then fit the distribution to a 2-dimensional quadratic form with 5 parameters de- 
scribing the expected elliptical shape: the coordinates (R.A. and Dec.) of the center of the 
ellipse, semi-major and -minor axis extents (a and /3), and the position angle of the ellips^. 
A 'quality' factor is evaluated to represent the goodness of the fit: it is the square root of the 
sum of the squares of the deviations for 8 points sampled along the contour where the value 
is expected to be 4.0, that is, 2 a from the maximum likelihood coordinates of the source. 

We quote the parameters of the ellipse that would contain 95% of the probability for the 
location of the source; for Gaussian errors this would be a radius of 2.45 a. An analysis of the 
deviations of 396 AGNs at high latitudes from the positions of the nearest LAT point sources 
indicated that the PSF width is underestimated, on average, by a factor of 1.10 ±0.05. Thus 
the final uncertainties reported by pointfit were scaled up by a factor of 1.1. To visually 
assess the fits, a Test Statistic map was made for each source, and these were considered in 
evaluating the analysis flags that are discussed in § 14. 8[ 

Twelve sources did not converge at the third step, converged to a point far away (> 1°) 
or were in crowded regions where the procedure (which does not have free parameters for the 
fluxes of nearby sources) may not be reliable. Those 12 were left at their gtfindsrc positions. 
They can be easily identified in the IFGL catalog because they have identical semimajor 
and semiminor axes for the source location uncertainty, and position angle 0. The LAT- 
detected pulsars and X-ray binaries, which were placed at the high-precision positions of 
these identified sources, have null values in the localization parameters. 

Figure [2] illustrates the resulting position errors as a function of the Test Statistic (TS) 
values obtained in § 14.31 The relatively large dispersion that is seen at a given TS is in part 
due to the local conditions (level of diffuse 7-ray emission) but primarily depends upon the 
source spectrum. Hard-spectrum sources are better localized than soft ones for the same TS 
(Fig. [3]) because the PSF is so much narrower at high energy. At our threshold of TS = 25 
the typical 95% position error is about 10', and most 95% errors are below 20'. 



■^In the FITS version of the IFGL catalog, a is CDNF_95_SEMIMAJ0R, f3 is CDNF_95_SEMIMIN0R, and (p is 
C0NF_95_PosAng; see Appendix D 
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Fig. 1. — Sky map of the LAT data for the time range analyzed in this paper, Aitoff 
projection in Galactic coordinates. The image shows 7-ray intensity for energies >300 MeV, 
in units of photons m~^ s^^ sr~^. 




Fig. 2. — 95% source location error (geometric mean of the two axes of the ellipse) as a 
function of Test Statistic (§ 14. 3p . The dashed line is a (TS)~^-^ trend for reference (not 
adjusted vertically). 
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Fig. 3. — 95% source location error multiplied by (TS)^'^ to remove the global trend (Fig. 
as a function of the photon spectral index from § 14.31 
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4.3. Significance and thresholding 

The detection and locahzation steps provide estimates of source significances. However, 
since the detection step does not use the energy information and the localization step fits 
only one source at a time, these estimates are not sufficiently accurate for use in the catalog. 
To better estimate the source significances we use a 3-dimensional maximum likelihood 
algorithm (gtlike) in unbinned mode (i.e., the position and energy of each event is considered 
individually) applied on the full energy range from 100 MeV to 100 GeV using the P6_V3 
IRFs (see § [2]). This is part of the standard Science Tools software package, currently at 
version 9rl5p5. The tool does not vary the source position, but does adjust the source 
spectrum. The underlying optimization engine is MinuitEI. The code works well with up to 
~30 free parameters, an important consideration for regions where sources are close enough 
together to partially overlap. The gtlike tool provides the best-fit parameters for each source 
and the Test Statistic TS = 2Alog(likelihood) between models with and without the source. 
The TS associated with each source is a measure of the source significance. Error estimates 
(and a full covariance matrix) are obtained from Minuit in the quadratic approximation 
around the best fit. For this stage we modeled the sources with simple power-law spectra. It 
should be noted that gtlike does not include the energy dispersion in the TS calculation (i.e., 
it assumes that the measured energy is the true energy). Given the 8 to 10% energy resolution 
of the LAT over the wide energy bands used in the present analyses, this approximation is 
justified. 

Because the fitted fluxes and spectra of the sources can be very sensitive to even slight 
errors in the spectral shape of the diffuse emission we allow the Galactic diffuse model (§ [3]) 
to be corrected (i.e., multiplied) locally by a power law in energy with free normalization 
and spectral slope. The slope varies between and 0.07 (making it harder) in the Galactic 
plane and the normalization by ± 10% (down from 0.15 and 20% for the BSL). The smaller 
excursions of that corrective slope when compared to the BSL reflect the better fit of the 
current diffuse model to the data. The normalization of the isotropic component of the 
diffuse emission (which represents the extragalactic and residual backgrounds) was left free. 
The three free parameters were separately adjusted in each Region of Interest (Rol). 

We split the sky into overlapping circular Rols. The parameters are free for sources 
in the central part of each Rol (Rol radius minus 7°), such that all free sources are well 
within the Rol even at low energy (7° is larger than rgs at 100 MeV). It is advantageous (for 
the global convergence over the entire sky) to use large Rols, but at the same time smaller 
Rols allow spectral variations of the diffuse emission relative to the model to be corrected 



'"^ http://lcgapp.cern.ch/project/cls/work-packages/mathlibs/minuit/doc/doc. html 
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in more detail. We set the Rol sizes so that not more than 8 sources are free at a time. 
Adding 3 parameters for the diffuse model, the total number of free parameters in each Rol 
is normally 19 at most. We needed 445 Rols to cover the 2433 seed positions. The Rol radii 
range between 9° and 15°. 

We proceed iteratively. All Rols are processed in parallel and a global current model 
is assembled after each step in which the best-fit parameters for each source are taken from 
the Rol whose center is closest to the source. The local model for each Rol includes sources 
up to 7° outside the Rol (which can contribute at low energy due to the broad PSF). Their 
parameters are fixed to their values in the global model at the previous step. The parameters 
of the sources inside the Rol but within 7° of the border are also fixed except in two cases 
(not considerered for the BSL analysis): 

• Sources within 2° of any source inside the central part, because they can influence the 
inner source. 2° is chosen to be larger than twice the containment radius at 1 GeV 
(2 X 0.8°) where the LAT sensitivity peaks (Fig. [T8|) . We leave both flux and spectral 
index free for these. 

• Very bright sources contributing more than 5% of the total counts in the Rol because 
they can influence the diffuse emission parameters. We leave only the flux free for 
these. 

All seed sources start at flux at the first step; the starting point for the slope is 2. We 
iterate over 5 steps; the fits change very little after the fourth. To facilitate the convergence 
the seed sources are not entered all at once. The brightest 10% of the sources are entered 
at the first step, 30% at the second step, and finally all at the third step. At each step we 
remove seed sources with low TS, raising the threshold for inclusion into the global model 
from 10 at the third step to 15 at the fourth and finally 25 at the last step. All seeds are 
reentered at the fourth step to avoid losing faint sources before the global model has fully 
converged. We have checked via simulations that removing the faint sources has little impact 
on the bright ones, much less than changing the diffuse model (§ 14. 6p . This procedure left 
1451 sources above threshold. The variation of the detection threshold across the sky and 
the dependence of the threshold on source spectrum are discussed in Appendix A. 

The TS of each source can be related to the probability that such an excess can be 
obtained from background fluctuations alone. The probability distributio n in such a situation 



TO 

(source over background) is not known precisely (jProtassov et al.ll2002l ). However since we 



consider only positive fluctuations, and each fit involves four degrees of freedom (two for 
position, plus flux and spectral index), the probability to get at least TS at a given position 



in the sky is close to 1/2 of the distribution with four degrees of freedom (IMattox et al. 
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1996i ). so that TS = 25 corresponds to a false detection probability of 2.5 x 10~^ or 4.1 a 
(one sided). For the BSL we considered only two degrees of freedom because the localization 
was based on a simpler algorithm which did not involve explicit minimization of the same 
likelihood function. 

The sources that we see are best (most strongly) detected around 1 GeV. This is approx- 
imately the median of the Pivot_Energy quantity in the catalog, i.e., the energy at which the 
uncertainties in normalization and spectral index for the power-law fit are uncorrelated. At 
1 GeV the 68% containment radius is approximately rgs = 0.8°. The number of independent 
elements in the sky (trials factor) is about 47r/(7rrgg) in which rgs is converted to radians. 
This is about 2 x 10^ so at a threshold of TS = 25 we expect less than 1 spurious source by 
chance only. If any, there might be a few very hard spurious sources in the catalog because 
hard sources have a smaller effective PSF so that the trials factor is larger. The main reason 
for potentially spurious sources, though, is our imperfect knowledge of the underlying diffuse 
emission (§ 14. 6p . 



4.4. Flux Determination 

The maximum likelihood method described in § l4.3l provides good estimates of the source 
significances and the overall spectral slope, but not very accurate estimates of the fluxes. 
This is because the spectra of most sources do not follow a single power law over that broad 
an energy range (three decades). Within the two most populous categories, the AGN often 
have broken power-law spectra and the pulsars have power-law spectra with an exponential 
cutoff. In both cases fitting a single power law over the entire range over-predicts the flux 
in the low-energy region of the spectrum, which contains the majority of the photons from 
the source, biasing the fluxes high. On the other hand the effect on the significance is low 
due to the broad PSF and high background at low energies. 

In addition, the significance is mostly obtained from GeV photons (Fig. [TS]) whereas 
the photon flux in the full range (above 100 MeV) is dominated by lower energy events so 
that the uncertainty on that flux can be quite large even for highly significant sources. For 
example, the typical relative uncertainty on the photon flux above 100 MeV is 23% for a 
TS = 100 source with spectral index 2.2. 

To provide better estimates of the source fluxes, we decided to spht the range into five 
energy bands from 100 to 300 MeV, 300 MeV to 1 GeV, 1 to 3 GeV, 3 to 10 GeV and 10 to 
100 GeV (the number of counts does not justify dividing the last decade into two bands). 
The list of sources remained the same in all bands. It is generally not possible to fit the 
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spectral index in each of those relatively narrow energy bands (and the flux estimate does 
not depend very much on the index), so we simply froze the spectral index of each source to 
the best flt over the full interval. The spectral bias to the Galactic diffuse emission (§ I4.3p 
was also frozen. 

The estimate from the sum of the flve bands is on average within 30% of the flux 
obtained from the global power-law flt (as described in § |2l with excursions up to a factor 
of 2. We have also compared those estimates with a more precise spectral model for the 
three bright pulsars (Vela, Geminga and the Crab). The sum of the flve fluxes is within 5% 
of the more precise flux estimate, whereas the power-law estimate is 25% too high for Vela 
and Geminga. However because it is not based on extrapolating a relatively well deflned 
power-law flt the relative uncertainty on that flux is even larger than that on the power-law 
flt, typically 50% for a TS = 100 source with spectral index 2.2. For that reason we do not 
show this very poorly measured quantity in Table [2l We provide instead the photon flux 
between 1 and 100 GeV (the sum of the three high energy bands), which is much better 
defined. The relative uncertainty on this flux is typically 18% for a TS = 100 source with 
spectral index 2.2. 

In contrast, the energy flux over the full band is better deflned than the photon flux 
because it does not depend as much on the poorly-measured low-energy fluxes. So we provide 
this quantity in Table [2j Here again the sum of the energy fluxes in the flve bands provides 
a more reliable estimate of the overall flux than the power-law flt. The relative uncertainty 
on the energy flux between 100 MeV and 100 GeV is typically 26% for a TS = 100 source 
with spectral index 2.2. 

An additional difficulty that does not exist when considering the full data is that, because 
we wish to provide the fluxes in all bands for all sources, we must handle the case of sources 
that are not signiflcant in one of the bands. Many sources have TS < 10 in one or several 
bands: 1135 in the 100 to 300 MeV band, 630 in the 300 MeV to 1 GeV band, 359 in the 
1 to 3 GeV band, 503 in the 3 to 10 GeV band and 800 in the 10 to 100 GeV band. There 
are even a number of sources which have upper limits in all bands, even though they are 
formally signiflcant (as deflned in § 14. 3p over the full energy range. It is particularly difficult 
to measure fluxes below 300 MeV because of the large source confusion and the modest 
effective area of the LAT at those energies with the current event cuts (§[2]). For the sources 
with poorly-measured fluxes (where TS < 10 or the nominal uncertainty of the flux is larger 
than half the flux itself), we replace the flux value from the likelihood analysis by a 2 a 
upper limit, indicating the upper limit by a in the flux uncertainty column of Tabled the 
corresponding columns of the FITS version of the IFGL catalog are described in Appendix 
D. The upper limit is obtained by looking for 2Alog(likelihood) = 4 when increasing the flux 
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observed at 8.4 GHz as part of VLA program AH996 ("Finding and Rejecting Associations for Fermi-hAT 7-ray sources"). This table is published in its entirety in the electronic edition of 
the Astrophysical Journal Supplements. A portion is shown here for guidance regarding its form and content. 
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from the maximum-likelihood value. When the maximum-likelihood value is very close to 
(i.e., the flux that maximizes the likelihood would be negative), solving 2Alog(likelihood) 
= 4 tends to underes timate t he up per limit. Whenever TS < 1 we switch to the Bayesian 



method proposed by iHelend (Il983l ). We do not use that method throughout because it is 



about five times slower to compute. 

The five fluxes provide a rough spectrum, allowing departures from a power law to be 
judged. This is the main advantage over the BSL scheme which involved only two bands. 
Examples of those rough spectra are given in Figures H] and for a bright pulsar (Vela) and 
a bright blazar (3C 454.3). In order to quantify departures from a power-law shape, we 
introduce a CurvatureJndex 

where i runs over all bands and Ff^ is the flux predicted in that band from the global power- 
law fit. fl^^ reflects the relative systematic uncertainty on effective area described in § 14.61 
It is set to 10, 5, 10, 15 and 20% in the bands [0.1,0.3], [0.3,1], [1,3], [3,10] and [10,100] GeV 
respectively. Note that this systematic uncertainty on the effective area is not included in 
the uncertainties reported in Table [3] (or in the FITS file), because this systematic factor 
cancels when comparing each of the band fluxes between different sources. We use for Fi 
and (Tj the best-fit and 1 a estimates even when the values are reported as upper limits in 
the table, both for computing the CurvatureJndex and the sums (photon flux and energy 
flux). 

Since the power-law fit involves two parameters (normalization and spectral index), C 
would be expected to follow a distribution with 5 — 2 = 3 degrees of freedom if the power- 
law hypothesis was true. At the 1% confidence level, the spectral shape is significantly 
different from a power law if (7 > 11.34. That condition is met by 225 sources (at 1% 
confidence, we expect 15 false positives). The curvature index is by no means an estimate of 
curvature itself, just a statistical indicator. A faint source with a strongly curved spectrum 
can have the same curvature index as a bright source with a slightly curved spectrum. Since 
the relative uncertainties on the fluxes in each band are quite different and depend on the 
spectral index itself, it is difficult to build a curvature indicator similar to the fractional 
variability for the light curves. The curvature index is also not exclusively an indicator of 
curvature. Any kind of deviation from the best fit power-law can trigger that index, although 
curvature is by far the most common. 
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Tabic 3. First LAT Catalog: Spectral Information 



Name IFGL 


r 


Ar 


Curv. 


100 MeV - 


300 MeV 


300 MeV - 


1 GeV 


1 GeV - 3 GeV 


3 GeV - 10 GeV 


10 GeV - 


100 GeV 














AFs'' . 






^Fa" - 




Fb<= AFs'^ 




J0000.8+6600C 


2.60 


0.09 




8.3 


0.0 


2.6 


1.9 


0.3 


6.4 


2.6 


0.6 


5.2 


3.7 


1.5 


4.1 


0.8 


0.0 


0.0 


JOOOO. 9-0745 


2.41 


0.20 




2.9 


0.0 


2.5 


0.5 


0.0 


3.0 


0.7 


0.0 


2.2 


3.8 


0.0 


3.7 


1.5 


0.0 


2.4 


JOOOl.9-4158 


1.92 


0.25 




2.1 


0.0 


1.8 


0.2 


0.0 


1.1 


0.6 


0.0 


2.2 


2.9 


1.1 


6.0 


1.7 


0.0 


0.0 


J0003. 1+6227 


2.53 


0.10 


T 


3.0 


0.0 


0.0 


1.7 


0.3 


6.8 


2.0 


0.5 


5.0 


4.3 


0.0 


1.4 


1.4 


0.0 


1.9 


J0004.3+2207 


2.35 


0.21 




1.8 


0.0 


0.9 


0.4 


0.0 


1.8 


0.4 


0.2 


3.2 


1.9 


0.9 


3.8 


0.9 


0.0 


0.0 


J0004.7-4737 


2.56 


0.17 




2.4 


0.8 


3.3 


0.3 


0.1 


3.6 


0.8 


0.2 


5.0 


2.0 


0.0 


2.0 


1.5 


0.0 


0.2 


J0005. 1+6829 


2.58 


0.12 




3.9 


0.0 


0.6 


1.3 


0.3 


5.3 


2.2 


0.0 


3.1 


4.9 


0.0 


2.0 


1.0 


0.0 


0.0 


J0005. 7+3815 


2.86 


0.13 




3.3 


0.9 


3.7 


0.5 


0.1 


4.4 


1.1 


0.0 


3.1 


2.5 


0.0 


1.6 


1.1 


0.0 


0.0 


J0006. 9+4652 


2.55 


0.11 




2.9 


0.9 


3.3 


0.7 


0.1 


6.4 


0.8 


0.3 


3.7 


3.0 


1.3 


4.3 


1.5 


0.0 


2.1 


J0007.0+7303 


1.97 


0.01 


T 


20.9 


1.2 


19.8 


12.0 


0.3 


60.5 


49.0 


1.3 


82.3 


135.6 


6.5 


57.6 


8.5 


1.6 


15.3 


J0008.3+1452 


2.00 


0.21 




0.9 


0.0 


0.0 


0.2 


0.0 


0.5 


0.6 


0.2 


3.8 


2.0 


0.9 


4.5 


1.2 


0.0 


0.0 


J0008.9+0635 


2.28 


0.22 




1.6 


0.0 


0.3 


0.5 


0.0 


3.1 


0.5 


0.0 


1.0 


2.2 


1.0 


4.0 


1.5 


0.0 


2.9 


J0009. 1+5031 


2.41 


0.13 




4.0 


0.0 


2.3 


0.5 


0.1 


4.6 


0.9 


0.3 


4.6 


3.5 


1.2 


5.1 


1.4 


0.0 


2.0 


JOOll. 1+0050 


2.51 


0.15 




1.2 


0.0 


0.0 


0.4 


0.1 


4.8 


0.5 


0.2 


4.1 


1.7 


0.0 


0.3 


0.9 


0.0 


0.0 


J0013.1-3952 


2.09 


0.22 




1.2 


0.0 


0.1 


0.3 


0.0 


1.4 


0.9 


0.0 


2.9 


2.3 


0.0 


1.4 


2.1 


0.0 


4.3 



Note. — This table is published in its entirety in the electronic edition of the Astrophysical Journal Supplements. A portion is shown 
here for guidance regarding its form and content. 

ajn units of 10~* photons cm~^ 

'°ln units of 10~® photons cm~^ s~^ 

"^In units of 10"^" photons cm~^ 
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Fig. 4. — Sample spectrum of Vela (IFGL J0835.3— 4510) generated from the five energy- 
band flux measurements in the catalog and plotted as EfAFi/ AEi, with Ei chosen to be 
the center of the energy bin in log space. The energy range of the integration is indicated 
by a horizontal bar. The vertical bar indicates the statistical error on the flux. The point 
at which these bars cross is not the same as the differential power per unit log bandwidth, 
E'^dF/dE at E^. The dashed lines (nearly coincident for this very bright source) reflect the 
uncertainties on the flux and index of the power-law fit to the full energy range in § 14.31 
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1 FGL J2253.9+1 608 - 3C 454.3 
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Fig. 5.— Spectrum of the bright blazar 3C 454.3 (IFGL J2253.9+1608). 
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4.5. Variability 

Variability is very common at 7-ray energies (particularly among accreting sources) and 
it is useful to estimate it. To that end we derive a variability index for each source by 
splitting the LAT data into a number of time intervals and deriving a flux for each source 
in each interval, using the same energy range as in § 14.31 (100 MeV to 100 GeV). We split 
the full 11-month interval into Vint = 11 intervals of about one month each (2624 ks or 
30.37 days). This is much more than the week used in the BSL, in order to preserve some 
statistical precision for the majority of faint sources we are dealing with here. It is also far 
enough from half the precession period of the orbit (~ 0.5 x 53.4 = 26.7 days) that we do not 
expect possible systematic effects as a function of off-axis angle to be coherent with those 
intervals. 

To avoid ending up with too large error bars in relatively short time intervals, we froze 
the spectral index of each source to the best fit over the full interval. Sources do vary in 
spectral shape as well as in flux, of course, but we do not aim at characterizing source 
variability here, just detecting it. It is very unlikely that a true variability in shape will be 
such that it will not show up in flux at all. In addition, little spectral vari ability was found 



in bright AGN where it would be detectable if present (lAbdo et al.ll2010kl ). Because we do 
not expect the diffuse emission to vary, we freeze the spectral adjustment of the Galactic 
diffuse component to the local (in the same Rol) best-fit index from the full interval. So the 
fitting procedure is the same as in § 14.41 with all spectral shape parameters frozen. 

The variability index is defined as a simple criterion: 

Fwt = (3) 
V=Y,'^i{F^-F^tf (4) 

i 

where i runs over the 11 intervals and cTj is the statistical uncertainty in Fj. As for the BSL we 
have added in quadrature a fraction f^^i = 3% of the flux for each interval Fj to the statistical 
error estimates ai (for each 1-month time interval) used to compute the variability index0. 
Since the weighted average flux F„t is not known a priori, V is expected, in the absence 
of variability, to follow a distribution with 10 (= Vint— 1) degrees of freedom. At the 



^In the FITS version of the IFGL catalog, the Flux_History and Unc_Flux_History columns contain Fi 
and af + {frciFiY, respectively; see Table [TT] 
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1% confidence level, the light curve is significantly different from a fiat one if \^ > 23.21. 
That condition is met by 241 sources (at 1% confidence, we expect 15 false positives). For 
those sources we provide directly in the FITS version of the table the maximum monthly 
fiux (Peak_Flux) and its uncertainty, as well as the time when it occurred (Time_Peak); see 
Table [11] for the column specifications. 

As in § 14.41 it often happens that a source is not significant in all intervals. To preserve 
the variability index (Eq. H]) we keep the best-fit value and its estimated error even when 
the source is not significant. This does not work, however, when the best fit is close to zero 
because in that case the log(likelihood) as a function of fiux is very asymmetric. Whenever 
TS < 10 or the nominal fiux uncertainty is larger than half the fiux itself we compute the 
2 a upper limit and replace the error estimate for that interval ((Xj) with half the difference 
between that upper limit and the best fit. This is an estimate of the error on the positive 
side only. Because the parabolic extrapolation often exceeds the log(likelihood) profile at 2 
a this is more conservative than computing the 1 a upper limit directly. The best fit itself 
is retained. Note that this error estimate can be a large overestimate of the error on the 
negative side, particularly in the deep Poisson regime at high energy. This explains why 
(Tj/Fj can be as high as 1 even when TS is 4 or so in that interval. As in § 14.41 we switch to 
the Bayesian method whenever TS < 1. 

Examples of light curves are given in Figures |6] and [7] for a bright constant source (the 
Vela pulsar) and a bright variable source (the blazar 3C 454.3). With the 3% systematic 
relative uncertainty no pulsar is found to be variable. The very brightest pulsars (Vela and 
Geminga) appear to have observed variability below 3%, so this may be overly conservative. 
It is not a critical parameter though, as it affects only the very brightest sources. 

The fractional variability of the sources is defined from the excess variance on top of 
the statistical and systematic fiuctuations: 



The typical fractional variability is 50%, with only a few strongly variable sou rces beyond 



SF/F = 1. This is qualitatively similar to what was reported on Figure 8 of lAbdo et al. 



( l2009nf ). The criterion we use is not sensitive to relative variations smaller than 60% at TS 
= 100. That limit goes down to 20% as TS increases to 1000. We are certainly missing many 
variable AGN below TS = 100 and up to TS = 1000. There is no indication that fainter 
sources are less variable than brighter ones; we simply cannot measure their variability. 



Both the curvature index and the variability index highlight certain types of sources. 
This is best illustrated on Figure |8] in which one is plotted against the other for the main 
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Fig. 6. — Light curve of Vela (IFGL J0835.3— 4510) for the 11-month interval analyzed 
for the IFGL catalog. The fluxes are integrated from 100 MeV to 100 GeV using single 
power-law fits and the error bars indicate the 1 o statistical errors. The grey band shows 
the time- averaged flux with the conservative 3% systematic error that we have adopted for 
evaluating the variabihty index. Vela is not seen to be variable even at the level of the 
statistical uncertainty. The spectrum of Vela is not well described by a power law and the 
fluxes shown here overestimate the true flux, but the overestimate does not depend on time. 
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1 FGL J2253.9+1 608 - 3C 454.3 
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Fig. 7. — Light curve of 3C 454.3 (IFGL J2253. 9+1608), which exhibits extreme variabihty. 
The grey band is the same 3% systematic uncertainty that we have adopted for evaluating 
the variabihty index. The triangles on the left and right indicate the value of the weighted 
average flux F^t that minimizes V in Eq. |H Owing to the systematic uncertainty term, for 
bright, highly- variable sources F„t can differ from the time- averaged flux (which we derive 
from a power-law fit to the integrated data set). 
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Fig. 8. — Variability index plotted as a function of curvature index (§ \4A\\ . The horizontal 
dashed line shows where we set the variable source limit, at > 23.21. The vertical 
dashed line shows where the spectra start deviating from a power-law, at C > 11.34. The 
cross standing out as very signifi cantly curved and variable is the source associated with 
LS I +61 303 kbdo et al.lboOQkh . 
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types of identified or associated sources (from the association procedure described in § |6]). 
One can clearly separate the pulsar branch at large curvature and small variability from the 
blazar branch at large variability and smaller curvature. 



4.6. Limitations and Systematic Uncertainties 

In this work we did not test for or account for source extension. All sources are assumed 
to be point-like. This is true for the major source populations in the GeV range (blazars, 
pulsars). On the other hand the TeV instruments have detected many extended sources 
in the Galactic plane, mostly pulsar wind nebulae and supernova remnants (SNRs), (e.g. 



Aharonian et al.l 120051) and the LAT has already started detecting extended sources (e.g. 
Abdo et al.l l2009il ). Because measuring extension over a PSF which varies so much with 
energy is delicate, we are not yet ready to address this matter systematically across all the 
sources in a large catalog such as this. 

We have addressed the issue of systematics for localization in § 14.21 Another related 
limitation is that of source confusion. This is of course strong in the inner Galaxy (§ 14. 7p 
but it is also a significant issue elsewhere. The average distance between sources outside 
the Galactic plane is 3° in IFGL, to be compared with a per photon containment radius 
res = 0-8° at 1 GeV where the sensitivity is best. The ratio between both numbers is not 
large enough that confusion can be neglected. The simplest way to quantify this is to look at 
the distribution of distances between each source and its nearest neighbor (-D„) in the area of 
the sky where the source density is approximately uniform, i.e., outside the Galactic plane. 
This is shown in Figure M The source concentration in the Galactic plane is very narrow 
(less than 1°) but we need to make sure that those sources do not get chosen as nearest 
neighbors so we select |6| > 10°. The histogram of Dn (after taking out the geometric factor 
as in Figure [9]) should follow 

H{Dn) = A^truePsrc exp (-VrD^^Psrc) (6) 

where psrc is the source density (number of sources per square degree) and A^true is the true 
number of sources (after correcting for missed sources due to confusion). The exponential 
term is the probability that no nearest source exists. It is apparent that, contrary to expec- 
tations, the histogram falls off toward Dn = 0. This indicates that confusion is important, 
even in the extragalactic sky. The effect disappears only at distances larger than 1.5°. To 
get Atrue, one may solve for the number of observed sources at distances beyond 1.5°. Since 
Psrc = Atrue/Aot in which Atot is the sky area at \b\ > 10°, this amounts to solving 



Arobs(> 1.5°) = A^true exp (-A^trueA/Aot) 



(7) 
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Fig. 9. — Distribution of the distances Dn to the nearest neighbors of all detected sources at 
|6| > 10°. The number of entries is divided by lirDn AD„ in which AZ}„ is the distance bin, in 
order to eliminate the 2-dimensional geometry. The overlaid curve is the expected Gaussian 
distribution for a uniform distribution of sources with no confusion (Eq. [6] normalized using 
Eq.ED. 
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in which Aq is the area up to 1.5°. This resuhs in A^true — ^obs = 80 missed sources on top 
of the A^obs = 1043 sources observed at |6| > 10°. Those missed sources are probably the 
reason for some of the asymmetries in the TS maps discussed in § 14.21 The conclusion is 
that globally we missed nearly 10% of the extragalactic sources. But because of the worse 
PSF at low energy soft sources are comparatively more affected than hard sources. This is 
approximately indicated by the difference between the full and the dashed lines on Figure [ 



Another important issue is the systematic uncertainties on the effective area of the 
instrument. At the time of the BSL we used pre-launch calibration and cautioned that 
there were indications that our effective area was reduced in flight due to p ile-up. Since 



then, the pile-up effect has been integrated in the simulation of the instrument (iRando et al. 



20091 ) and many tests have shown that the resulting calibration (P6_V3) is consistent with 
the data. The estimate of the remaining systematic uncertainty is 10% at 100 MeV, 5% 
at 500 MeV rising to 20% at 10 GeV and above. This uncertainty applies uniformly to 
all sources. Our relative errors (comparing one source to another or the same source as a 
function of time) are much smaller, as indicated in § 14.51 The fluxes resulting from this new 
calibration are systematically higher than the BSL fluxes. For example, the fluxes of the 
three brightest pulsars (Vela, Geminga and Crab) are about 30% larger in IFGL than in 
the BSL. The differences are more pronounced for soft sources than hard ones. This implies 
also that the IFGL flu xes are significantly larger than the EGRET fluxes in the 3EG catalog 
flHartman et al.lll999h which happene d to be close to the BSL fluxes. As shown by diffuse 



(lAbdo et al.ll2009hl ) and point source ( lAbdo et al.ll2009ii l2010el ) observations, the LAT data 



produce spectra systematically steeper than those reported in EGRET analysis. LAT fluxes 
are greater at energies below 200 MeV and less at energies above a few GeV. 

The model of diffuse emission is the other important source of uncertainties. Contrary 
to the effective area, it does not affect all sources equally: its effects are smaller outside the 
Galactic plane {\b\ > 10°) where the diffuse emission is faint and varying on large angular 
scales. It is also less of a problem in the high energy bands (> 3 GeV) where the PSF 
is sharp enough that the sources dominate the background under the PSF. But it is a 
serious issue inside the Galactic plane (|6| < 10°) in the low energy bands (< 1 GeV) and 
particularly inside the Galactic ridge (|/| < 60°) where the diffuse emission is strongest and 
very structured, following the molecular cloud distribution. It is not easy to assess precisely 
how large the uncertainty is, for lack of a proper reference model. We discuss the Galactic 
ridge more specifically in § 14.71 For an automatic assessment we have tried re-extracting 
the source fluxes assuming a different diffuse model, derived from GALPROP (as we did 
for the BSL) but with protons and electrons adjusted to the data (globally). The model 
reference is 54_87Xexph7S. The results show that the systematic uncertainty more or less 
follows the statistical one (i.e., it is larger for fainter sources in relative terms) and is of the 
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same order. More precisely, the dispersion is 0.7 o" on flux and 0.5 a on spectral index at 
\h\ > 10°, and 1.8 cr on flux and 1.2 a on spectral index at |6| < 10°. We have not increased 
the errors accordingly, though, because this alternative model does not fit the data as well 
as the reference model. From that point of view we may expect this estimate to be an upper 
limit. On the other hand both models rely on nearly the same set of H I and CO maps of the 
gas in the interstellar medium, which we know are an imperfect representation of the mass. 
That is, potentially large systematic uncertainties are not accounted for by the comparison. 
So we present the figures as qualitative estimates. 

4.7. Sources Toward Local Interstellar Clouds and the Galactic Ridge 

Figure [10] shows an example of the striking, and physically unlikely, correspondence 
between the IFGL sources and tracers of the column density of interstellar gas, in this case 
E(B-V) reddening. The sources in Orion appear to be tightly associated with the regions 
with greatest column densities. Yet no particular classes of 7-ray emitters are known to be 
associated with interstellar cloud complexes. Young SNRs would be resolved in the radio and 
in 7-rays in the nearby clouds outside the Galactic plane. Even if radio-quiet pulsars were 
the sources, they would not be expected to be aligned so closely with the regions of highest 
column densities. The implication is that peak column densities are being systematically 
underestimated in the model for the Galactic diffuse emission used in the analysis. E(B-V) 
is not directly used in the model; as described in § |3]the column densities are derived from 
surveys of H I and CO line emission, the latter as a tracer of molecular hydrogen. An E(B- 
V) 'residual' map E(B-V)i.es; representing interstellar reddening that is not correlated with 
A^(H I) or W{CO), is included in the model. So the peak column densities would need to 
be underestimated both in CO and E(B-V). We are studying the effect and strategies for 
validating the model for Galactic diffuse emission at high column densities. 

In addition to the concerns about the accuracy of column densities toward the peaks of 
interstellar clouds, self absorption of H I at low latitudes can introduce small angular-scale 
underestimates of the column densities, and intensities of the diffuse emission. The current, 
half-degree binned, model of the interstellar emission used for the source analysis also cannot 
capture structure on smaller angular scales. Bright structure on finer scales could be detected 
as unresolved point sources. 

Figure [To] illustrates another, better understood issue with the model for Galactic diffuse 
emission. In the Orion nebula (near /, 6 ~ 209°, —19.5°), a massive star-forming region that 
is extremely bright in the infrared, the infrared color corrections used to evaluate E(B- 
V) are inaccurate and the column densities inferred from E(B-V) are underestimated; the 
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depression in E(B-V) in the nebula does not correspond to a decreased column density 
of gas. The model of Galactic diffuse emission was constructed with E(B-V)i.cs allowed 
to be a signed correction for the column densities inferred from H I and CO lines. This 
made essential improvements in large regions, but in discrete directions toward massive star- 
forming regions, negative E(B-V)res can introduce deep depressions in the predicted diffuse 
emission. Figure [11] illustrates the depression around the S225 star-forming region, and its 
close correspondence with a IFGL source. 

Similar considerations relate to the sources at low latitudes in the inner Galaxy. The 
density of unassociated sources in the Galactic ridge (300° < / < 60°, |6| < 1°) is very 
high (Fig. [T2|) . and their latitude distribution is exceedingly narrow (Fig. [T3|) . If these IFGL 
sources are true 7-ray emitters they must have a very small scale height in the Milky Way, like 
that of the yo ungest massive star- forming regions, traced by ultracompact H II regions (~25' 



FWHM, e.g., iGiveon et al.ll2005l ). or be quite distant and hence very luminous. The IFGL 
sources do not have an obvious correspondence with the ultracompact H II regions, and the 
latter are not plausible 7-ray sources, but owing to the effects described above, embedded 
star-forming regions can influence tracers of gas and dust and thereby potentially introduce 
small-scale errors in the model of Galactic diffuse emission. The inferred luminosities of the 
IFGL sources in the Galactic ridge would be quite high if the scale height of their distribution 
is characteristic of most tracers of Population I objects. For a relatively narrow dispersion of 
40 pc, the characteristic distances of these sources are ~11 kpc (i.e., more distant than the 
Galactic center) and the 7-ray luminosities exce ed 10^^ erg s~^ (i. e., more than an order of 



magnitude more luminous than the Vela pulsar, lAbdo et al.ll2009il ). For broader dispersions 



about the plane, the distances and luminosities would increase correspondingly. 

The IFGL sources toward the peaks of local interstellar clouds and the Galactic ridge 
all have analysis flags set (§ 14. 8p in the catalog. We have also added a designator 'c' to 
their names to indicate that they are to be considered as potentially confused 
with interstellar diffuse emission or perhaps spurious. In addition, the 'c' designator 
is used for unidentified IFGL sources in crowded regions of high source density outside the 
Galactic ridge, as a caution about the complications due to PSF overlaps The 'c' desig- 
nator, thus applied to 161 of the IFGL sources, is a warning that the existence 
of the source or its measured properties (location, flux, spectrum) may not be 
reliable. 
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Fig. 10. — Overlay of IFGL sources on a square root color scale representation of E(B- 
V) reddening in Orion ( iSchlegel et al.lll998[ l. The units for the color bar are magnitudes 
of reddening. The white circles indicate the positions and 95% confidence regions for the 
IFGL sources in the field. The magenta circles indicate the effective (spectrally-weighted) 
68% containments for photons >500 MeV associated with each source that is positionally 
correlated with the clouds; these circles can be considered to represent the region of the 
sky most relevant for the definition of each source. The yellow contours around (/, b) = 
(-151°, -19.2°), (-153.5°, -16.2°), and (-146.3°, -12.8°) outline the regions with negative 
reddening residuals caused by errors in the dust infrared color corrections near young clusters 
of IR sources. The Orion Nebula is near L 6 ~ 209°, —19.5°. 
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Fig. 11. — Overlay of IFGL J0541. 1+3542 sources on a square root color scale representa- 
tion of E(B-V)res residual reddening in the S235 H II region. The units for the color bar 
are magnitudes of reddening. The small circle indicates the position and 95% confidence re- 
gion for the IFGL source. The large circle indicates the effective (spectrally-weighted) 68% 
containment for p hotons with energies >500 MeV. The W{CO) intensity contours (from 
Dame et al.l 120011 ) at 6, 15, and 30 K km s~^ trace the column density of molecular gas. 
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Fig. 12. — Overlay of IFGL sources (magenta circles) on an image of the intensity W{CO) 
of the 2.6 mm line of CO (jPame et al.ll200ll ). in a segment of the Galactic ridge in the 
first quadrant. Th e black asterisks indicate the positions of ultracompact H II regions 
( iGiveon et al.ll2005l ). which are similarly narrowly distributed about the Galactic equator. 
For the IFGL sources, the inner circles indicate the 95% confidence regions for the locations 
and the outer circles the approximate extents of the spectrally-weighted PSF for energies 
>1 GeV. The latter give an approximate sense of the 'regions of infiuence' for the IFGL 
sources. 
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Fig. 13. — Latitude distribution of unassociated/unidentified IFGL sources in the Galactic 
ridge (300° < / < 60°). 
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4.8. Analysis Flags 

We have identified a number of conditions tliat can slied doubt on a source. Tliey are 
described in Table HI As noted, setting of flags 4 and 5 depends on the energy band in which 
a source is detected. The high energy bands are always more favorable for source detection 
and characterization, as a result of the narrower PSF at high energy, so the flags are set on 
the basis of the highest band in which a source is significant. Flag 5 signals confusion and 
depends on a reference distance 6'ref- Because statistics are better at low energy (enough 
events to sample the core of the PSF), ^ref is set to the Full Width at Half Maximum 
(FWHM) there (minimum distance to have two peaks with a local minimum in between in 
the counts map). At high energy there are fewer events so 6'ref is set to the larger value 2rQs. 
In the intermediate bands we interpolate between FWHM and 2r68- In the FITS version of 
the catalog, these flags are summarized in a single integer column (Flags; see Appendix D). 
Each condition is indicated by one bit among the 16 bits forming Flags. The bit is raised 
(set to 1) in the dubious case, so that good sources have Flags = 0. 

5. The IFGL Catalog 

In this section we tabulate the quantities listed in Table [2] for each source; see Table [5] 
for descriptions of the columns. The source designation is IFGL JHHMM.m+DDMM where the 1 
refers to this being the first LAT catalog, FGL represents Fermi Gamma-ray LAT. Sources 
close to the Galactic ridge and some nearby interstellar cloud complexes are assigned names 
of the form IFGL JHHMM.m+DDMMc, where the c indicates that caution should be used in 
interpreting or analyzing these sources. Errors in the model of interstellar diffuse emission, 
or an unusually high density of sources, are suspected to affect the measured properties or 
even existence of these sources (see § 14.71) . 

The designations of the classes that we use to categorize the IFGL sources are listed in 
Table [6] along with the numbers of sources assigned to each class. We distinguish between 
associated and identified sources, with associations depending primarily on close positional 
correspondence (see § 16. 2p and identifications requiring measurement of correlated variability 
at other wavelengths or characterization of the IFGL source by its angular extent (see § l6.3p . 
Sources associated with SNRs are often also associated with PWNs and pulsars, and the 
SNRs themselves are often not point-like. We do not attempt to distinguish among the 
possible classifications and instead in Table [7] list plausible associations of each class for 
unidentified IFGL sources that are found to be associated with SNRs. 

The summed photon flux for 1-100 GeV (-F35; the subscript ij indicates the energy range 



Table 4. Definitions of tlie Analysis Flags 



Flag^ Meaning 

1 Source with TS > 35 which went to TS < 25 when changing 
the diffuse model (§ 14. 6p . Note that sources with TS < 35 
are not flagged with this bit because normal statistical 
fluctuations can push them to TS < 25. 

2 Moved beyond its 95% error ellipse when changing the diffuse 
model. 

3 Flux or spectral index changed by more than 3 a when 
changing the diffuse model. Requires also that the flux change 
by more than 35% (to not flag strong sources). 

4 Source-to-background ratio less than 30% in highest band in 
which TS > 25. Background is integrated over vrrgg 

or 1 square degree, whichever is smaller. 

5 Closer than ^ref from a brighter neighbor, ^^ef is defined in the 
highest-energy band in which source TS > 25. 6'ref is set to 
2.6° (FWHM) below 300 MeV, 1.52° between 300 MeV and 

1 GeV, 0.84° between 1 GeV and 3 GeV, and 2r68 above 3 GeV. 

6 On top of an interstellar gas clump or small-scale defect in the 
model of diffuse emission. 

7 Unstable position determination; result from gtfindsrc outside 
the 95% ellipse from pointlike (see § 14. 2p . 

8 pointlike did not converge. Position from gtfindsrc. 

9 Elliptical quality > 10 in pointlike (i.e., TS contour does not 
look elliptical). 



^In the FITS version the values are encoded in a single column, with Flag 
n having value 2*^"~^^ For information about the FITS version of the table 
see Appendix D and § [51 
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as 10* - 10^' MeV) and the energy flux for 100 MeV to 100 GeV in Table [2] are evaluated 
from the fluxes in bands presented in Table [31 The energy fluxes in each band are derived 
on the assumption that the spectral shape is the same as for the overall power-law fit (r25). 
This is an approximation but the bands are narrow enough that the energy fluxes are not 
strongly dependent on the spectral index. We do not present the integrated photon flux 
for 100 MeV to 100 GeV because the relative uncertainties tend to be very large, because 
of the limited effective area in the lower energy bands. Figure [H] shows the distribution of 
integrated fluxes (1 GeV - 100 GeV) for all of the sources in the catalog. The dynamic range 
is approximately 2.5 orders of magnitude, owing both to the intrinsic range of fluxes of the 
7-ray sources as well as their spectral hardnesses. 

The procedure used to designate associated sources and class designations is described 
in § 16.21 Figure US] shows the distribution of the IFGL sources on the sky by source class, 
and Figure [16] shows just the inner Galaxy. 

Figure [T71 selects just the variable sources, i.e., those with <1% chance of being steady 
sources in the monthly light curves and those with spectra incompatible with power laws. 
The variable sources are seen to be predominantly outside the Galactic plane, and many 
are associated with blazars. The spectrally-curved sources have a distribution much more 
confined to the Galactic equator. 

The electronic version of the IFGL catalog, available in FITS format from the Fermi 
Science Support Center, includes all of the information in these tables plus the monthly light 
curves from which the variability index values were derived and pivot energies for the overall 
power-law fits; see Table [Til 

6. Source Association and Identification 

Even with the good angular resolution of LAT, source location accuracy is typically 
not precise enough to make a firm identification based on positional coincidence alone. A 
typical LAT error region contains numerous stars, galaxies. X-ray sources, infrared sources, 
and radio sources. Determination of the nature of a given LAT source must therefore rely 
on more information than simply location, including time variability, spectral information, 
and availability of sufficient energy and a plausible physical process at the source to produce 
7-rays. 

In this analysis, we make a clear distinction between a source identification and an 
association with an object at another wavelength. A firm identification of a source is based 
on a timing characteristic such as a periodicity for a pulsar or binary or a variability correlated 



Table 5. LAT First Catalog Description 



Column Description 



Name IFGL JHHMM.m+DDMM[c] , constructed according to lAU Specifications for Nomenclature; m is decimal 
minutes of R.A.; in the name R.A. and Decl. are truncated at 0.1 decimal minutes and 1', respectively; 
c indicates that based on the region of the sky the source is considered to be potentially confused 
with Galactic diffuse emission 

R.A. Right Ascension, J2000, deg, 3 decimal places 

Decl. Declination, J2000, deg, 3 decimal places 

/ Galactic Longitude, deg, 3 decimal places 

b Galactic Latitude, deg, 3 decimal places 

a Semimajor radius of 95% confidence region, deg, 3 decimal places 

/9 Semiminor radius of 95% confidence region, deg, 3 decimal places 

(p Position angle of 95% confidence region, deg. East of North, decimal places 

o" Significance derived from likelihood Test Statistic for 100 MeV— 100 GeV analysis, 1 decimal place 

F35 Photon flux for 1 GeV-100 GeV, 10"^ ph cm s , summed over 3 bands, 1 decimal place 

AF35 l-CT uncertainty on i^35 , same units and precision 

S25 Energy flux for 100 MeV-100 GeV, 10~^^ erg cm~-^ s~^, from power-law fit, 1 decimal place 

AS25 1-f uncertainty on 525, same units and precision 

r Photon number power-law index, 100 MeV-100 GeV, 2 decimal places 

AT 1 a uncertainty of photon number power-law index, 100 MeV— 100 GeV, 2 decimal places 

Curv. T indicates < 1% chance that the power-law spectrum is a good fit to the 5-band fluxes; see note in text 

Var. T indicates < 1% chance of being a steady source; see note in text 

Flag See Table 1 for deflnitions of the fiag numbers 

7-ray Assoc. Positional associations with OFGL, 3EG, EGR, or AGILE sources 

TeV Positional association with a TeVCat source, P for angular size <20', E for extended 

Class Like 'ID' in 3EG catalog, but with more detail (see Table (6)1. Capital letters indicate firm identifications; 

lower-case letters indicate associations. 

ID or Assoc. Designator of identified or associated source 

Ref. Reference to associated paper (s) 
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Table 6. LAT IFGL Source Classes 



Description 


Designator 


Number Assoc. (ID) 


Pulsar, X-ray or radio, identified by pulsations 


psr (PSR) 


7 (56) 


Pulsar, radio quiet (LAT PSR, subset of above) 


PSR 


24 


Pulsar wind nebula 


pwn (PWN) 


2 (3) 


Supernova remnant 


t (SNR) 


41 (3) 


Globular Cluster 


glc (GLC) 


8(0) 


Micro-quasar object: X-ray binary (black hole 


mqo (MQO) 


0(1) 


or neutron star) with radio jet 






Other X-ray binary 


hxb (HXB) 


0(2) 


BL Lac type of blazar 


bzb (BZB) 


295 (0) 


FSRQ type of blazar 


bzq (BZQ) 


274 (4) 


Non-blazar active galaxy 


agn (AGN) 


28 (0) 


Active galaxy of uncertain type 


agu (AGU) 


92 (0) 


Normal galaxy 


gal (GAL) 


6(0) 


Starburst galaxy 


sbg (SBG) 


2 (0) 


Unassociated 




630 



Note. — The designation 'f indicates potential association with SNR or PWN (see 
Table 0. Designations shown in capital letters are firm identifications; lower case 
letters indi cate associations. In the case of AGN, many of the associations have high 
confidence llAbdo et al.ll2010lh . Among the pulsars, those with names beginning with 
LAT were discovered with the LAT. For the normal galaxy class, 5 of the associations 
are with the Large Magellanic Cloud. In the FITS version of the IFGL catalog, the f 
designator is replaced with 'spp'; see Appendix D. 
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Table 7. Potential Associations for Sources Near SNRs 



Name IFGL 



Assoc. 



JUzl0.5+DiOOC 


133 


.01 





,82 


SNR 


/"I -1 on '7 1 r\-\ 


.3 






JUizU.U+OZO ( 


132 


.80 


1 


,80 


blNH 


-i Qo T 1 m 
Uloz. ( "hUl, 


.3 






JUoUU.l+52d7 


155. 


.63 


6. 


.34 


SNR 


Cj15d.2+05. 


.7 






JU5Uo.z+45zo 


161 


.64 


2. 


.35 


SNR 


(jrlbU.9+Uz, 


.6 






JUOOO.O+Z i ii 


180 


.59 


—2 


.16 


blNR 


1 on n m 


,7 






JUooo.y+oiUo 


179 


.08 


2 


.65 


olNK 


1 vn Ci 1 no 


,6 






J Udod.UH-U45oc 


206 


.74 


— 1 


.15 


olNK 


(^onc; 1 nn 


,5 






Tnooo o ,10,10 

JUozo.o— 4248 


260. 


.37 


—3, 


.15 


SNR 


(j2d0.4— 03. 


,4 






JUo41.y— 4ozU 


265. 


.18 


—2. 


.58 


bJNK 


Cjr^OO.y— Uo. 


.3 






JUcS04.U— 40oz 


2DD. 


.64 


— 1. 


.09 


oJNK 


LtZOD.z — Ui. 


.2 


(Vela Junior) 




JiUlo.D — OOOO 


2o4. 


.32 


— 1 


.70 


oJNK 


/~1r)Oy1 Q ni 

tjrZo4.o— Ui. 


.8 


(MSH 10-53) 




J iiiy.4 — OiZ ( C 


292. 


.1 / 


— 


.52 


OTVT'D 

blNR 


i^ono o nn 
CjjZyz.z— UU. 


,5, 


PWN G292.2- 


0.5; 


Ti 1 Q/1 Q i^nc; c; 
J iio4.e5— dUoo 


onQ 

zyo. 


■7 "7 
. ( / 


U. 


,0 t 


GM"D 
olNK. 


r~^onQ c 1 nn 


,D, 


PWN G293.8+0.6 


Jlzlo.7— Dz4Uc 


298. 


.63 


— 


,12 


SNR 


/^nno r' nn 

Cj298.6— 00. 


,0 






Tio^o T o o n „ 

Jld4d.7— Dzciyc 


308. 


.89 


— 0, 


.39 


SNR 


ono o nn 

(j308.8— 00. 


,1 






Tiorrr\ o 

JIOOU.O — Dzl^C 


309. 


.80 


— 


.12 


SNR 


onn o 1 nn 

CjrciUy.o+UU. 


.0 






J i4iU.O— DlzoC 


312. 


.21 


— 0. 


.03 


oJNK 


Ooiz.4— UU. 


■4, 


PSR J1410-6132 


Ji442.U— DUUUc 


316. 


OA 

.34 


— u 


.05 


bJNK 


/~i QIC Q nn 
Lrolb.o— UU. 


.0 






JloUi.D— 4zU4 


327. 


.30 


14, 


,55 


SNR 


00*7 1 14 


,6 






T-iri/i '7 c;ni'7 

J ioi4. /— oyi ( 


320. 


.33 


— 1 


,35 


olNR 


/^Qon /I ni 


,2 


(MSH 15-52) 




J iOZi.5— o4c 


322. 


.03 


— 


,38 


blNR 


/^Qoi n nn 


,3 






Jlo52.4— 5oU9 


326. 


.25 


— 1 


.71 


SNR 


o o ni 

(j326.3— 01. 


.8, 


PWN G326.3- 


1.8 


JlolY.o— oluoc 


332. 


.39 


— 


.40 


SNR 


/"looo yi nn 

(jr332.4— 00 


.4 






J1640.8-4634C 


3.38, 


,29 


-0 


.06 


SNR 


G338.3-00. 


,0, 


PWN G338.3- 


0.0 


J1711.7-3944C 


347, 


,15 


-0 


.19 


SNR 


G347.3-00. 


,5 






J1725. 5-2832 


357, 


,96 


3. 


,91 


SNR 


G358.0+03. 


,8 






J1745.6-2900C 


359. 


.94 


-0, 


.05 


SNR 


GOOO.0+00. 


,0, 


PWN G359.95 


-0.1 


J1801.3-2322C 


6. 


.57 


-0. 


.22 


SNR 


G006.4-00. 


.1 


(W28) 




J1805.2-2137C 


8. 


.55 


-0. 


.14 


SNR 


G008.7-00. 


.1 


(W30) 




J1806.8-2109C 


9. 


.13 


-0 


.24 


SNR 


G008.7-00. 


.1 


(W30) 




J1834.3-0842C 


23. 


,24 


-0, 


.20 


SNR 


G023.3-00. 


,3 


(W41) 




J1834.7-0709C 


24. 


,67 


0. 


,43 


SNR 


G024.7+00. 


,6 






J1916.0+1110C 


45. 


,67 


-0, 


.31 


SNR 


G045.7-00. 


,4 






J2020.0+4049 


78. 


.37 


2. 


.53 


SNR 


G078.2+02. 


.1 






J2042.3+5041 


88. 


.66 


5. 


.19 


SNR 


G089.0+04. 


.7 






J2046.0+4954 


88. 


.42 


4. 


.24 


SNR 


G089.0+04. 


.7 






J2046.4+3041 


73. 


.40 


-7, 


.79 


SNR 


G074.0-08. 


.5 






J2049. 1+3142 


74. 


,56 


-7, 


.60 


SNR 


G074.0-08. 


,5 






J2055. 2+3144 


75, 


,40 


-8, 


.59 


SNR 


G074.0-08. 


,5 






J2057.4+3057 


75, 


,11 


-9, 


.46 


SNR 


G074.0-08. 


,5 






J2323.4+5849 


111. 


,74 


-2, 


.12 


SNR 


Gill. 7-02. 


.1 







Note. — See text. These sources are marked with a f in Table 6. They may be pulsars 
or PWN rather than the SNR named. 




Fig. 14. — Distribution of the fluxes of the IFGL sources in the energy range 1 GeV - 
100 GeV. 
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□ No association ^ Possible association witli nearby SNR or PWN 
X AGN - blazar * Starburst Gal * Pulsar * Pulsar w/PWN 
X AGN - unknown + Galaxy <;> PWN a Globular cluster 

X AGN - non blazar o SNR k XRB or MQO 



Fig. 15. — The 1451 IFGL catalog sources, showing locations on the sky (in Galactic coor- 
dinates with Aitoff projection) and associated source class, coded according to the legend. 
Gray symbols indicate unassociated sources, blue indicate associated sources, and red sym- 
bols firmly identified sources (primarily pulsars). For this plot the bzb and bzq designators 
have been combined ( "AGN-blazar" ) , as have hxb and mqo ("XRB or MQO"). The sources 
possibly associated with SNR, PSR or PWN (those indicated by a dagger in Table [2]) are 
listed as "Potential SNR". 
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Fig. 16. — The IFGL catalog sources in the inner Galactic region. Sources are indicated in 
the same manner as Figure [151 
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with observations at another wavelength in the case of a blazar, or on measurement of finite 
angular extent, which is the case for some Galactic sources, e.g., SNRs. An association is 
defined as a positional coincidence that is statistically unlikely to have occurred by chance 
between a plausible 7-r ay-producing object and a LAT source. 

For the IFGL catalog, the approach to designating associations and identifications in- 
volves three steps: 

1. A test for statistically-si gnificant detectio n s of c lasse s of sources, b ased on a defined 
protocol as described by iTorres fc Reimerl (120051 ) and iReimerl (120071 ) has been carried 
out. Based on these results, some potential source classes are deemed unlikely even if 
individual sources in such classes might be positionally coincident with IFGL sources. 

2. A gen eral automated sou rce association analysis, enhanced from the version used for the 
BSL (lAbdo et al.ll2009nl ). has been applied to the sources in the catalog. This method 
relies principally, although not exclusively, on comparing the local space density of 
plausible source classes with the number of positional associations found for a class 
with IFGL sources. 

3. For pulsars and binary systems with clearly-identified periodic emission, a firm identi- 
fication can be established. 



Each of these methods is described below. 



6.1. Protocol for population identification 



The idea for a protocol for population testi ng in 7-ray source data s uch as represented 
by the Fermi Ist-year catalog was introduced by lTorres &: Reimerl (120051 ). a paper to which 
we refer for general background. The test was devised to provide high levels of confidence in 
population classifica t ion w ith small number statistics. Essential to this protocol is the use of 
Feldman &: CousinsI (119981 ) confidence level intervals in a priori, physically selected samples 
of plausible 7-ray emitters. 

Thus, we consider spatial correlations between the Ist-year catalog sources (1451 detec- 
tions) and a priori selected (on physical grounds or earlier hints in 7-ray data of previous 
missions) sets of astrophysical sources, details on which are given below. To test for spatial 
correlations, each of the Fermi sources is described by a centroid position and an uncertainty. 
The latter is described by an ellipse, with semi-minor and semi-major axis and a position 
angle; with all these values being taken from the catalog at 95% CL. Those sources that are 
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identified beyond doubt by timing (all 7-ray pulsars and three X-ray binaries) have their 
positions assigned to be coincident with the cor responding astrop hysical sources (as located 



in other wavelengths or by LAT pulsar timing; ISmith et al.ll2008l ). However, for population 
searches, we maintain the uncertainty of their detections around these positions. In order 
to test systematics we also analyze the case in which all IFGL sources have their positional 
uncertainties enlarged by increasing their corresponding semi-minor and semi-major axes by 
20%. This is done consistently both for the real set of IFGL sources, and for each of the 
simulated sets generated by Monte Carlo (see below), and is meant as a check of the results 
obtained with nominal uncertainty values. 

Let C{A) represent the number of coincidences between candidate counterparts for pop- 
ulation A and LAT sources. In case of finding the same Fermi detections spatially coincident 
with several astrophysical objects of the same class, we count all of the coincidences as one. 
An example of this is seen for millisecond pulsars pertaining to the same globular cluster, 
e.g., 47 Tuc. To determine the number of excess coincidences above the noise level, i.e., the 
sources that are expected to correlate by chance, we subtract the background b caused by 
these random coincidences, given the number and distribution of LAT sources and members 
of the testing population. 

To obtain this latter number, we produce Monte Carlo simulations shuffling the orien- 
tation of the elliptical position uncertainty (but not its size) and the centroid position of 
each of the actual IFGL sources, thus generating sets of fake LAT source detections. Each 
simulated IFGL catalog is constructed to have the total number of sources (1451 detections), 
and both the longitude and latitude distributions of the real set. In order to obtain the most 
conservative results we proceed to simulate fake LAT source catalogs by maintaining the 
latitude (longitude) histogram of the real set with 5, and 10 (30, 60) degrees binning, and 
subsequently take the largest value for the expected average number of random positional 
coincidences among them (and this is called b). 

The number of excess coincidences above chance associations is then S{A) = C{A)~b{A). 
This number is used to test the null hypothesis: Population A is not 'j-ray emitting at flux- 
levels detectable by the Ist-year LAT catalog. The predicted number of source coincidences 
for this hypothesis is equal to 0, and the total expected events if this hypothesis is valid is 
equal to b. The greater the excess of the real number of correlations over the corresponding 
fe-value for that population, the easier it is to rule out the null hypothesis. 

The testing power of a sample of finite size is limited: if using the same set of data, 
claiming the discovery of one population affects the level of confidence by which one can 
claim the discovery of a second. In order to control the reliability of our results, we require 
that the combination of all of our claims be bounded by a probability of 10~^, which then 
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becomes the total budget B. This low probability provides an overall significance of about 
5cr, which implies individual claims of populations must arise with higher confidence. The 
total budget can then be divided into individual ones for each population, P4, Pg, etc., such 
that ^ ■ Pi = B. Then, population A will be claimed as detected with x% CL if and only if: 



the Poissonian cumulative probability (CP) for obtaining the real number of spatial 
coincidences (or more) as a result of chance coincidence is less than the a priori assigned 
budget Pa (as opposed to being less than only the larger, total budget) and 

the number of excesses £ {A) is beyond the upper lirn it of the corresponding confidence 
interval for x% CL, from lFeldman fc Cousins! (119981 ). 



The latte r value, x%, for eac h popu lation is then the confidence level obtained using 
the tables in iFeldman fc CousinsI ( 1l998l ) for which the upper end of the interval equals 
the real number of excesses. (With number of events 6, background 6, recall that in the 
null hypothesis, there are expected events above background.) The values of Pi for each 
population are chosen very conservatively: 0.01% B to each of the only two populations 
that were unambiguously identified in the EGRET catalog, pulsars and blazars; 0.1% B 
to millisecond pulsars and EGRET-coincident SNRs, which were hinted at in the EGRET 
catalog; and the rest of the budget equally distributed into those classes where no high-energy 
7-ray emission was previously reported. 

Finally, we provide brief notes on the populations selected for the test. The selections 
were made before the start of Fermi operations, as described in the internal LAT-document 
AM-090670, which lists each specific source selected in each class. 



Blaza rs: The selected blazars are a subset of the CGRABS catalog ( jHealey et al. 



20081 ). which is a unifor m all-sky survey of EGRET-like blazars, selected by their 



Figure-of-merit (FOM: see lHealey et al.ll2008l ). In total the CGRABS catalog includes 
1625 sources. To assemble the blazar test hst we cut the CGRABS catalog at the 
small est FoM value for w hich all high-confidence EGRET blazars of the 3EG cat- 
alog ( Hartman et al. [l999[), which are also Mattox-blazars ( Mattox et al. 2001), are 



included (see ISowards-Emmerd et al.l l2003l . for a comparison of 3EG and "Mattox" 



blazars). This cut is at FoM=0.111, and results in a total of 215 sources and constitutes 
the blazar list used for the application of the protocol. 



^Available from https://oraweb.slac.stanford.edu/pls/slacquery/DOCUMENTS.STARTUP?PROJECT=GLAST 



- 53 - 



Misaligned jet sour ces: We start with the compilation of extragalactic jet sources 
(iLiu &: Zhangll2002l ) which has 661 entries collected from the literature (as of December 
2000). This list contains radio galaxies, radio quasars, BL Lac objects and Seyfert 
galaxies. For a redshift cut at z < 0.032 (chosen to select a moderate number of the 
closest members of this class) the list features 51 sources: 34 are clasified as radio 
galaxies (including all 7-ray radio galaxies detected with high confidence so far), 16 as 
Seyfert galaxies, and one as BL Lac object (not yet detected in 7-rays and not included 
in the CGRABS list). 3C 111 and 3C 120 are added on the basis of evidence for their 
detection by EGRET. 



Starbu rsts: This list includes all starburst galaxies for which in the study of lTorres et al. 
(I2OO4J ) the combination of gas content, cosmic-ray density, and distance indicated a flux 
detectable by the LAT in one year. 

Halo Dwarf Galaxies: The list consists of presently-known objects in this category up 
to 1000 kly. This sample of Milky Way satellite galaxies is useful for probing localized 
7-ray excesses due to annihilation of dark matter particles. Dwarf spheroidal galaxies 
(dSph) may be manifestations of the largest clumps predicted by the CDM scenario. 

Galaxy clusters: We restrict this sample to 30 clusters with the highest values of mass- 
to- dis t ance-sq uared , M/(P, from the HIFLUCS sample of the iReiprich fc Bohringer 
( I2OO2I ) sample, complemented by clusters that are interesting individually (Bullet clus- 
ter, RX J1347.5), as well as those reported to have radio haloes. Galaxy clusters with 
reported hard X-ray emission, and those predicted from large scale struc ture formation 
simulations to be detectable by the LAT in one year (jPfrommeii l2008l ) are implicitly 
included the sample. 



Pulsars: These are selected from the ATNF catalog (IManchester et al.l 120051 ). with a 



cut ow E > 10^^ erg s ^ . 

Millisecond Pulsars: These are also selected from the ATNF catalog to have periods 
less than 10 ms, with the same cut as above in E. 

Magnetars: We include all 13 known SGRs and AXPs. 

EGRET-SNRs: These are all S NRs that were found to be spatially coincident with an 
EGRET source, as discussed by lTorres et al.l (l2003al ). 



TeV shell-type SNRs: these were intended to be tested separately to answer the ques- 
tion whether these objects also emit GeV 7-rays on at least the level of the LAT 1 
yr sensitivity. The sample comprised 4 sources (RX J1713.7— 3946, RX J0852— 460, 
RCW 86 and Gas A). 
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Star-star binaries: These are from the Vllth WR catalog of Ivan der Huchtl (l200ll ). 
requiring 0.001 Ekin^tot/ {.'^T^d\) > 10, resulting in a sample of 41 sources. The factor 
0.001 is motivated by a reasonable 10% acceleration efficiency, i.e., wind-energy-to- 
relativistic-particle-energy-conversion efficiency, and 1% radiative efficiency in these 
environments, while the factor 10 relates to the LAT sensitivity anticipated at the 
time of protocol application. 

Star-compact object binaries: These are all known 7-ray binaries and microquasars 
(MQs). 

Binary pulsars: We include all objects known in this category. 

Globular clusters: Globular clusters are known to contain a relatively large number of 
millisecond pulsars (MSPs) whose individual and collective emission in the X-ray and 7- 
ray energy bands may be detectable by the LAT. Here the aim is to test for the collective 
emission from MSPs, since MSPs are principally able to accelerate leptons at the shock 
waves originated in collisions of the pulsar winds and/or inside pulsar magnetospheres, 
and inside a globular cluster these are subsequently able to Comptonize stellar and 
microwave background radiation. The globular clusters in this list are restricted to be 
closer than 6 kpc. 



6.1.1. Protocol results 



Table |H] summarizes the results of the application of the protocol for population iden- 
tification to the IFGL catalog. The first column shows the name of the test population. 
The second column {Np) shows the number of astrophysical objects included in the test of 
the population. The third column shows the expected number of random coincidences (or 
background, b{A)) between the 1451 LAT sources in the IFGL catalog, and the Np candi- 
dates in each population as explained above. The fourth column (C) gives the actual number 
of positional coincidences between the members of each population and the IFGL catalog. 
The fifth column shows the cumulative probability (CP) of obtaining C or more positional 
coincidences between IFGL sources and the astrophysical candidates purely by chance. The 
sixth column shows the a priori assigned probability {P{A)) for discovery of each popula- 
tion, the sum of all P giving the total budget for population discovery. The seventh column 
answers yes or no to the question of whether the actual random probability for an equal 
or a larger number of excesses to occur by chance (CP) is larger than the a priori assigned 
budget for this population. For the cases in which this answer is yes, the last column gives 
the confidence level of the detection of the corresponding population; alternatively stated. 
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n No association n Possible association with nearby SNR or PWN 
X AGN - blazar * Starburst Gal * Pulsar * Pulsar w/PWN 
X AGN - unknown + Galaxy ^ PWN a Globular cluster 

X AGN - non blazar o SNR a XRB or MQO 
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-90 



n No association n Possible association with nearby SNR or PWN 
X AGN - blazar * Starburst Gal * Pulsar * Pulsar w/PWN 
K AGN - unknown + Galaxy <=> PWN a Globular cluster 

X AGN - non blazar o SNR a XRB or MQO 



Fig. 17. — The IFGL catalog sources with variable (left) and spectrally-curved {right) 
sources highlighted. Sources with the variability flag or spectral curvature flags set are 
shown according to their associated class, as red symbols, coded as marked as in the legend 
while sources showing no evidence of variability or with spectra compatible with power laws 
are shown as black dots. 



Table 8. Results from Application of the Population Protocol 



Test Population Np b C CP P (CP < P)? CL 



Galactic Populations 



Pulsars 


215 


1.440 


30 


5.3 


X 


10" 


-29 


1 





X 


IQ- 


-9 


yes 


> 99.999% 


Millisecond Pulsars 


23 


0.050 


7 


1.5 


X 


10- 


-13 


1 





X 


10- 


-8 


yes 


99.89% 


EGRET SNRs 


23 


1.590 


13 


1.5 


X 


10- 


-8 


1 





X 


10- 


-8 


no 




TeV SNRs 


4 


0.920 


3 


6.6 


X 


10- 


-2 


9 


9 


X 


10- 


-7 


no 




Magnetars 


13 


0.120 













9 


9 


X 


10- 


-7 






WR- binaries 


41 


0.260 













9 


9 


X 


10- 


-7 






MQ/7-ray bin. 


17 


0.140 


3 


4.1 


X 


10- 


-4 


9 


9 


X 


10- 


-7 


no 




Binary pulsars 


10 


0.040 













9 


9 


X 


10- 


-7 






Globular clusters 


29 


0.240 


4 


1.1 


X 


10- 


-4 


9 


9 


X 


10- 


-7 


no 




jxtragalactic Populat 


ions 




























Blazars 


215 


0.480 


61 


0.0 








1 





X 


10- 


-9 


yes 


> 99.999% 


Misaligned jet sources 


53 


0.150 


5 


5.5 


X 


10- 


-7 


9 


9 


X 


10- 


-7 


yes 


99.25% 


Starbursts 


15 


0.050 


4 


2.5 


X 


10- 


-7 


9 


9 


X 


10- 


-7 


yes 


97.89% 


Galaxy clusters 


48 


0.150 













9 


9 


X 


10- 


-7 






Dwarf spheriodals 


18 


0.070 













9 


9 


X 


10- 


-7 
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this is the confidence with which we are ruhng out the hypothesis that such a population is 
not present among the LAT detections (the null hypothesis). 

Not surprisingly, the source populations already conclusively identified in the EGRET 
era are found with the highest confidence in the investigated IFGL coincidences even when 
very strict thresholds were chosen for associations to be claimed, justifying a posteriori the 
very low budget assigned to these populations. 

To study the sensitivity of the results to the sizes of the error ellipses for the sources, 
we also evaluated the coincidences and chance probabilities when the extents of the ellipses 
are increased by a (very conservative) 20%. The findings for misaligned jet sources and 
EGRET SNR populations were not robust against this change. When the sizes of the source 
location regions are increased the probabilities of chance associations necessarily increase 
for any given source. For the misaligned jet source population, this was enough to push 
the population below the a priori probability budget. For the EGRET SNRs, enlarging the 
source location regions had the offsetting effect of adding one more coincidence with IFGL 
sources so the net result was a chance probability below the a priori threshold. 

Spatial coincidences found between LAT sources and individual millisecond PSRs (MSPs) 
and those tested via their globular cluster environment are related in the following sense. 
Whereas for an isolated MSP there are not many alternative scenarios for producing de- 
tectable 7-ray emissio n, and in fact sever al of them have been detected and identified by 



their 7-ray per iodicitv (lAbdo et al.ll2009b[) . in the case of coincidences with globular clusters 
(e.g., 47 Tuc, lAbdo et al.ll2009dl ). an ambiguity exists as long as no pulsed emission from 
one of MSPs in the globular cluster is found. Thus, the aim was to search not only for the 
existence of a population of MSPs, but also to distinguish between environments in which 
MSPs are found as populations as well. The fact that the class of globular clusters is not 
detected as a population but the millisecond pulsars are can be interpreted as related to 
the size of the population of MSPs in a globular cluster. Based on the membership of the 
IFGL catalog apparently only those clusters hosting a large number of MSPs were found 
coincident with catalog sources (e.g., Ter 5, 47 Tuc, M28), too small a number to claim them 
as a population of 7-ray sources in the framework of this test. 



The previously-repor ted relation between EGRET sources and SNRs (e.g., ISturner fc Dermer 



I995I : [Torres et al.l l2003bl ) cannot be confirmed with IFGL sources using the present test. 
The chance probability of the number of coincidences seen with IFGL sources is very small, 
but the a priori probability budget assigned for this population was smaller still (Table |H]). 
We note also that pulsars have been detected by the LAT in some of the SNRs in our sample. 



The fact that the misaligned jet sources are found as a population supports the individual 
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IFGL source association findings, although the test did not distinguish between radio galaxies 
and Seyferts, and does not show a significant correlation with ranking via radio core flux. 
It is to be noted that this class is found only when the uncertainties in the positions of the 
LAT sources are not enlarged, when the cumulative probability is low, but found to be just 
below the a priori assigned budget for the class. 

Starburst galaxies arc detected as a population in LAT data, in numbers that add up to 
the noted individual detections of M82 and NGC 253, and suggests the potential for future 
discoveries. 

Particularly interesting is the non-detection of galaxy clusters, which is not only a 
statistically significant result, but even a case of zero coincidences. With the absence of 
even a single galaxy cluster coincidence from the tested sample, we can conclude that X-ray 
bright nearby galaxy clusters and those exhibiting a radio halo do not constitute a source 
population above the IFGL sensitivity limit: the models that we used to select the candidate 
galaxy clusters overestimate the energy conversation into particle acceleration at GeV and 
greater energies. 

Dwarf spheroidals are not found coincident with LAT sources either. The result is com- 
patible with the null hypothesis that this population does not exist above the sensitivity 
limit of the catalog, although in some dark-matter based models of 7-ray production they 
would have been expected to be seen at the sensitivity level of the IFGL catalog. The Galac- 
tic populations of magnetars, binary pulsars, and WR binaries are similarly not detected. 
The latter is particularly interesting; whereas positional associations of WR stars with LAT 
sources can be found using the extensive list of known WR stars, those having the largest 
wind energies do not present any correlation and thus we disregard WRs as 7-ray emitters 
at the sensitivity limit of the IFGL catalog. 

Finally, we note that whereas we found three sources correlated with MQ-7-ray binaries 
(all of them secured by timing) the population as such is not detected. We can conclude 
that the binaries identified in the IFGL catalog are certainly special objects, but probably 
not archetypal for a population of similar objects. Similarly, the same occurs for TeV SNRs: 
we find some correlated objects, but not enough to claim a population discovery due to the 
large expected background (b) for this population. 

6.2. Automated Source Associations 

Our approach for automated source association follows closely that used for the BSL, 
although we enlarged our database of catalogs of potential counterparts and improved our 
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calibration scheme to control more precisely the expected number of fa lse associations. The 



association procedure, which follows essentially the ideas developed by iMattox et al.l (119971 ) 
for the identification of EGRET sources with flat-spectrum radio sources, is described in 
detail in Appendices B and C. Here we only summarize the essential steps of the automated 
source association procedure. 

The automated source association is based on a list of catalogs that contain potential 
counterparts of LAT sources. This list has been compiled based either on prior knowledge 
about classes of high-energy 7-ray emitters or on theoretical expectations. In total, 32 
catalogs (some of which are subselections from 24 primary catalogs, see Table [9]) have been 
searched for counterparts covering AGNs (and in particular blazars), nearby and starburst 
galaxies, pulsars and their nebulae, massive stars and star clusters, and X-ray binaries. For 
the BSL analysis, only 14 catalogs were searched. In addition, some of the catalogs have 
been enlarged somewhat (e.g., BZCAT, CRATES, and SNR total; cf. Table [9] with Table 1 
of the BSL). Furthermore, since millisecond pulsars have now been established as sources 
of 7-rays, we split the ATNF catalog into normal pulsars and millisecond pulsars. For the 
latter we require the pulse period P < 0.1 s and the period derivative P < 10~^^ s/s. This 
results in 139 objects. The remaining 1545 objects are considered as normal pulsars. We 
divided those into high and low E/(P categories, the latter being defined as E/d'^ < 5 x 10^^ 
erg kpc~^ s~^. This results in 84 high E/d"^ pulsars and 1461 low E/d"^ pulsars. Note 
that this energy-fiux selecting is unrelated to the selection on spin-down luminosity applied 
for the population study in § 16.11 For SNRs, we divided the Green catalog into two lists, 
one containing all objects that can be considered as point-like for the LAT (157 objects) 
and one containing extended supernova remnants (117 objects, diameters greater than 20'); 
these subsets are denoted Small and Large, respectively. We search for counterparts at 
radio frequencies using the VLBA Calibrator Survey and at TeV energies using the TeVCat 
catalog. For the latter we divide the TeVCat catalog into Small and Large angular size 
subsets at 40'. We also search for coincidences between IFGL sources and BSL (OFGL), 
AGILE, and EGRET sources. The complete list of catalogs, the numbers of objects they 
contain, and the references are presented in Table [91 
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Table 9. Catalogs Used for the Automatic Source Association and Results 







p 

1 




V 

- ^ ass 




\-Mals(-, 


\ 

/ 




LAT pulsars 


56 




0.1 


56 


n.a. 


0. 


,4 


1 


High E/(f pulsars 


84 





.024 


24 


n.a. 


0. 


,6 


2 


Low E/dP pulsars 


1461 





.011 


1 


n.a. 


0. 


,3 


2 


MiUisecond pulsars 


139 


0, 


.278 


20 


n.a. 


1. 


,0 


2 


Pulsar wind nebulae 


69 





.049 


27 


0.3 


0. 


,9 


1 


High-mass X-ray binaries 


114 





.010 


3 


n.a. 


0. 


,3 


3 


Low-mass X-ray binaries 


187 





.050 


8 


0.4 


0. 


,5 


4 


Small (<20') SNRs 


157 





.021 


11 


0.7 


0, 


,7 


5 


stars 


378 





.015 


1 


< 0.1 


< 0. 


,1 


6 


WR stars 


226 





.013 


11 


0.3 


0. 


,2 


7 


LBV stars 


35 





.026 


2 


0.3 


0. 


,6 


8 


Open clusters 


1689 


0, 


.013 


1 


0.1 


0. 


,4 


9 


Globular clusters 


147 


0, 


.272 


8 


< 0.1 


0. 


,5 


10 


Nearby galaxies 


276 





.066 


5 


0.4 


0. 


,4 


11 


Starburst galaxies 


14 




0.5 


2 


< 0.1 


< 0, 


,1 


12 


Blazars (BZCAT) 


2837 





.308 


487 


8.9 


6. 


,8 


13 


Blazars (CGRaBS) 


1625 





.238 


282 


4.7 


4. 


,1 


14 


Blazars (CRATES) 


11499 





.333 


490 


17.2 


17. 


,8 


15 


BL Lac 


1122 





.224 


218 


2.8 


2. 


,8 


16 


AGN 


21727 





.021 


11 


0.7 


0. 


,8 


16 


QSO 


85221 





.166 


147 


7.3 


4. 


,9 


16 


Seyfert galaxies 


16343 





.041 


24 


2.0 


1. 


,6 


16 


Radio-loud Scyfcrt galaxies 


29 




0.1 


4 


< 0.1 


< 0. 


,1 


1 


VLBA Calibrator Survey 


4558 





.266 


484 


11.5 


10. 


,1 


17 


Small (<40') TeV sources 


92 





.037 


42 


0.6 


0. 


,8 


18 


Large (>40') TeV sources^ 


11 




n.a. 


13 


n.a. 


7, 


,5 


18 


Large (>20') SNRs'f 


117 




n.a. 


48 


n.a. 


18. 


,1 


5 


Dwarf galaxies^ 


14 




n.a. 


7 


n.a. 


2. 


,1 


1 


1st AGILE catalog* 


47 




n.a. 


52 


n.a. 


18. 


,6 


19 


3rd EGRET catalog* 


271 




n.a. 


107 


n.a. 


25. 


,4 


20 


EGR catalog* 


189 




n.a. 


66 


n.a. 


9. 


,1 


21 
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For each catalog in the hst, we make use of Bayes' theorem to compute the posterior 
probabihties Pik that an object i from the catalog is the correct association of the LAT 
source k: 



ik 



\ -'prior / 



-^^^^rior ijllG prior probability that counterpart i is detectable by the LAT, Ofc and are the 
axes of the ellipse at 1 a, smaller than the semimajor and semiminor axes at 95% confidence 
by a factor ^—2 log(0.05) = 2.45, pk is the local counterpart density around source k, and 



/ coiP'id) — (hk) sin^(0 — 



for a given position angle (j) between LAT source k and the counterpart i, (pk being the 
position angle of the error ellipse, and r being the angular separation between LAT source k 
and counterpart i (see Appendix B). For each catalog, prior probabilities Pprior are assigned 
so that 

A^false= Yl (1-^^-^) (10) 

Pik^Pthr 

gives the expected number of false associations that have posterior probabilities above the 
threshold Pthr (see Appendix C). The corresponding prior probabilities are quoted in Table [9] 
(column 3) for all catalogs. 

For the automated association of the IFGL catalog we set Pthr = 0.8, which means that 
each individual association has a < 20% chance of being spurious. This is different from 
the approach we took for the BSL paper where we constrained the expected number of false 
associations for each catalog to A^faisc ^ 1? which imposed a relatively tight constraint on 
source classes with large numbers of associations (such as blazars and pulsars) while source 
classes with only few associations had a relatively loose constraint. Now, each individual 
association stands on an equal footing by ha ying a well defined probability for being spurious. 



We note that the First LAT AGN Catalog (lAbdo et al.l 1201011 ) applies the same association 
method as we use here but includes associations for AGNs down to Pthr = 0.1. 

For a number of catalogs in our list the Bayesian method cannot be applied since either 
(1) the location uncertainty of the counterpart is larger than the location uncertainty for 
the IFGL source, or (2) the counterpart is an extended source; see notes to Table M In 
the first case, we consider all objects i as associations for which the separation to the LAT 
source k is less than the quadratic sum of the 95% confidence error radius of counterpart i 
and the semimajor axis a^. For the bright EGRET pulsars Crab, Geminga, Vela, and PSR 
J1709— 4429, we also list associations, even though in many cases the EGRET-measured 
locations are formally inconsistent with the positions of the pulsars; there is no doubt that 
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Table 9 — Continued 



Name Objects Ppnor A^ass A^faise {N^se) Ref. 



Bright Source List (OFGL) 205 n.a. 195^ n.a. 3.9 22 



References. 



1 Collaboration internal; 2 Manchester et al 



Liu et al. 
J2004h: 



van der Hucht 



Massaro et al 



(1200 ih: 8 IClark et all fcOOSl): 9 



IClHarrisI fll996h: llSchmidt et al.l p998h: 1 2 IXhompson et al 



fl2009h : l^lHealev et all J2007h : 15 iHealev et al. 



(B); 



f|2006l): 4 Ihiu et all fcoO Th : 5 ICreenl toO^: 6 iMafz-Apellaniz et al. 



Dias et al. 



Veron-Cettv fc VeronI J2006h : 17 iKovalevI ([2009^ and ref. 



120021 ): 



( 2007 ): 13 
koo4 \: 16 
therein, 



http://astrogeo.Org/vlbi/solutions/2009c_astro/i 18 http: / /tevcat. uchicago.edu/ 
('Defa ult' and 'Newly Announ ced' categories, version 3. 1 00); 19 iPittori et al 



f 2009h: 20 lHartman et al.l (|l999l ): 2l lCasandjian fc Grenieil (|2008h : 22 lAbdo et al 
fl2009nh 



* Catalog for which the location uncertainties of the counterparts are greater 
than the location uncertainties for the IFGL sources. 

"''Catalog for which the counterparts are spatially extended sources. 

'^See § 16.2.101 
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EGRET detected these pulsars. In the second case, we assume that the counterparts have a 
circular extension and consider all objects i as associations for which the circular extension 
overlaps with a circle of radius around the LAT source k. 



The results of the automated association procedure for each of the external catalogs are 
summarized in Table [91 For each catalog, we give the number 



of LAT sources that have been associated with objects in a given catalog (column 4). Fur- 
thermore, we compute the expected number of false associations A^faise using Eq. flTU]) for 
those catalogs which have been associated with the Bayesian method (column 5). We can- 
not give meaningful results for pulsars and high-mass X-ray binaries since, for identified 
objects, the positions have been fixed in the catalog to their high- precision locations (§ H]), 
and consequently, their posterior probabilities are by definition 1. However, we alternatively 
estimated the expected number of false associations using Monte Carlo simulations of 100 
realizations of fake LAT catalogs, for which no physical associations with counterpart catalog 
objects are expected (see Appendix C). We quote the resulting estimates (iVfaise) in column 6 
of Table |9l We find A^faise ~ ( A^faise) which confirms that the posterior probabilities computed 
by the automatic association procedure are accurate (otherwise Eq. f ITU]) would not hold). 

In total we find that 821 of the 1451 sources in the IFGL catalog (56%) have been 
associated with a least one non-7-ray counterpart by the automated procedure at the 80% 
confidence level. 779 IFGL sources (54%) have been associated using the Bayesian method 
while the remaining 42 sources are spatial coincidences based on overlap of the error regions 
or source extents. From simulations we expect that 57.3 among the 821 sources (7%) are as- 
sociated spuriously. Considering only the Bayesian associations, 37.5 among the 779 sources 
(5%) are expected to be spurious. In the following we discuss the automated association 
results in some detail. Associations with TeV sources are discussed in § [71 



6.2.1. Automated association summary 




1 



(11) 



6.2.2. Blazars 



Our association procedure contains 4 catalogs to cover the blazar source class (BZCAT, 
CGRaBS, CRATES, BL Lac) and these catalogs have a substantial number of objects in 
common. In total we find 689 IFGL sources associated with sources from at least one of 
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the 4 blazar catalogs. 2 of these sources (IFGL J0047.3-2512 and IFGL J0956.5+6938) 
are the starburst galaxies M 82 and NGC 253 (both found in the CRATES catalog), and 
2 sources (IFGL J0319.7+4130 and IFGL J1325.6-4300) are the radio galaxies NGC 1275 
and Cen A. This leaves 685 blazar candidates among the IFGL sources. 

We further note that 282 of the IFGL sources associated with blazars also have counter- 
parts in the VL BA calibrator sur vey (VCS) which we added to our list of catalogs following 
the suggestion of lKovalevI (j2009al ) who found 111 associations for this catalog among the BSL 
sources. For 37 of the 484 IFGL sources associated with VCS objects, the VCS association 
is the only counterpart found among all catalogs. Most of these 37 sources are located at 
low Galactic latitudes, a region in which our 4 blazar catalogs are incomplete. Many of the 
low-lat itude VCS asso ciations thus may be related to blazars situated close to the Galactic 
plane ( iKovalevI l2009al ) . 



6.2.3. Other AGNs 



We find 24 IFGL sources that are associated with objects from the two Seyfert galaxy 
catalogs. Among those, only two sources are not also associated with blazars: IFGL J0840. 8+1310 
(3C 207.0) and IFGL J1230.8+1223 (M 87: lAbdo et aDl2009gh 3C 207.0 is a lobe-dominated 
quasar and and M 87 is a radio galaxy. 

147 of the IFGL sourc es are associated with AGNs and QSOs from the catalog of 
Veron-Cetty fc Veroru (120061 ). yet all of these are also associated with blazars. Therefore 
apparently most of the IFGL non-blazar associations are either sources that are in fact 
blazars, yet are not classified as such in our catalogs, or they are nearby radio galaxies. 
In particular, we do not find convincing evidence for coincidences of IFGL sources with 
non-blazar Seyfert galaxies. 



6.2.4. Normal Galaxies 



We find 2 associations with nearby starburst galaxies (lAbdo et al.ll2010al l IFGL J0047.3- 
2512 (NGC 253) and IFGL J0956 5+693 8 (M 82). Both galaxies have also been detected 
at TeV energies ( lAcero et al.ll2009al ) and ( IVERITAS Collaboration! 120091 ) and hence can be 
considered as high-confidence IFGL associations. 

Seven IFGL sources are found to coincide with dwarf galaxies: 6 are associated with the 
Large Magellanic Cloud (LMC), 1 is associated with the Small Magellanic Cloud (SMC), both 
galaxies being extended. These sources probably correspond to local maxima of extended 
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emission features and probably do not represent real point sources in the field. Regarding 
the LMC, in Table |2] only 5 of the sources are indicated as being associated with the LMC 
The exception is IFGL J0600.7— 7037 which is also associated with PKS 0601—70, a blaza r 
that was significantly variable during the time span of the IFGL catalog flAbdo et al.ll2009j ). 

Thus NGC 253, M 82, the LMC, and the SMC are so far the only normal galaxies that 
have been associated with sources in the IFGL catalog. We note that IFGL J1305.4— 4928 
is associated with NGC 4945, which is a starburst galaxy that is also classified as a Seyfert 
II AGN. 



6.2.5. Pulsars, pulsar wind nebulae and globular clusters 

56 IFGL sources have been identified as pulsars through their 7-ray pulsations. For 
these sources we list only the pulsar identification in Table HI This is not to be taken to 
mean that we have necessarily ruled out contributions from known or unknown PWNs or 
SNRs; in fact for the sources identified with the Crab and Vela pulsars we have also identified 
PWNs; these sources have both class assignments 

In addition to the 56 seen pulsating, we find 3 more associations with the high E/d"^ 
subset of pulsars from the ATNF catalog: 

• IFGL J1119.4-6127C (PSR J1119-6127) 

• IFGL J1410.3-6128C (PSR J1410-6132) 

• IFGL J1648.4-4609C (PSR J1648-4611) 



These 3 sources are good candidates for young energetic 7-ray pulsars, although IFGL J1119.4— 6127 
and IFGL J1410.3-6128c are also associated with SNRs and/or PWNs (see Table [7]). 

Among the IFGL sources that are associated with pulsar wind nebulae (PWNs), only 
6 are not also associated with known pulsars: 

• IFGL J1134.8-6055 (PWN G293.8+0.6) 

• IFGL J1552.4-5609 (PWN G326.3-1.8) 

• IFGL J1635.7-4715C (PWN G337.2+0.1) 

• IFGL J1640.8-4634C (PWN G338.3-0.0) 
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• IFGL J1745.6-2900C, the Galactic center source (PWN G359.95-0.04) 

• IFGL J1746.4-2849C (PWN GO. 13-0.11) 

It remains to be shown whether the LAT indeed detects these PWNs, or whether the 7-ray 
emission arises from the yet-unknown pulsars that power the nebulae, or potentially from 
an associated SNR. Because of the ambiguity we list only two positional associations with 
PWNs in Table El for IFGL J1635.7-4715c and IFGL J1746.4-2849c, which do not also 
have associations with known SNRs. The others are included in Table [71 

Among the 20 IFGL sources that are associated with millisecond pulsars, 11 are not 
associated with known 7-ray pulsars. Among those 11,5 are associated with globular clusters 
and the remaining 6 may indeed be Galactic field 7-ray millisecond pulsars: 

• IFGL J0610. 7-2059 (PSR J0610-2100) 

• IFGL J1024.6-0718 (PSR J1024-0719) 

• IFGL J1600.7-3055 (PSR J1600-3053) 

• IFGL J1713.9+0750 (PSR J1713+0747) 

• IFGL J1811.3-1959C (PSR J1810-2005) 

• IFGL J1959.6+2047 (PSR B1957+20) 

Finally, we find that 8 IFGL sources are associated with globular clusters. None of those 
have alternative associations different from millisecond pulsars or low-mass X-ray binaries 
(which both are known source populations residing in globular clusters), which makes the 
reality of these associations even more plausible. 



6.2.6. Supernova remnants 

Our automated association procedure associates 59 IFGL sources with SNRs. Not 
counting associations that also include pulsars detected by the LAT, the total is 41 (Tab. [7]). 
Of those, 5 are associated with small angular size (diameter less than 20') SNRs: 

• IFGL Jl 134.8-6055 (SNR G293.8+00.6, also associated with PWN G293.8+0.6) 

• IFGL J1213.7-6240 (SNR G298.6-00.0) 
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IFGL J1617.5-5105 (SNR G332.4-00.4) 

IFGL J1640.8-4634 (SNR G338.3-00.0, also associated with PWN G338.3-0.0) 



IFGL J2323.4+5849 (SNR Gill. 7-02.1, aka Gas A; lAbdo et al.ll2ninfl ) 



Except for IFGL J2323.4+5849 and possibly IFGL J1213.7-6240, the presence of alternative 
associations to PWNs or a low-mass X-ray binary make the physical association of these 
sources to SNRs questionable. 

Some of the associations with SNRs have already been suggested based on morphology 
analyses of the LAT sources: 



IFGL J1856.1+0122 (G034.7-00.4, aka W44: lAbdo et al.ll2010hl ) 



IFGL J1922.9+1411 (G049.2-00.7, aka W51G: lAbdo et al.ll2009ih 



IFGL J0617.2+2233 (G189.1+03.0, aka IG 443: lAbdo et al.ll2010il ) 



We consider these three SNRs to be identified sources in the IFGL catalo g; see ^16.31 Further 
interesting associations due to the presence of OH masers in the SNR (IHewitt et al.l 120091 ) 
are IFGL J1805.2-2137c and IFGL J1806.8-2109c (both overlapping G008.7-00.1, also 
known as W30). 



6.2.1. Association of IFGL J1745.6—2900c with the Galactic center 

With a position {l,b) = (359.941°, -0.051°) and a 95% confinement radius of 1.1', 
IFGL J1745.6— 2900c is the source closest to the Galactic center. In this direction, many 
catalogs contain objects and consequently we find a large number of formal associations to 
this source. Specifically, IFGL J1745.6— 2900c is formally associated with the pulsar wind 
nebula G359.95-0.04, the SNR GOOO.0+00.0, the VGS object J1745-2900, 4 low-mass X- 
ray binaries, 6 LBV stars, and 10 Wolf-Rayet stars and 2 TeV sources. We are unable to 
distinguish on the basis of our association scheme among these possibilities, although some 
are more plausible physically. Eventually, the spectral energy distribution of the source or 
any characteristic time-variability may help to narrow down the possibilities. 
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6.2.8. X-ray binaries 



Three IFGL sources have been identified by their orbital modulations as high-mass 
X-ray binaries (HMXB): 



IFGL J0240.5+6113 (LS I+6r303: lAbdo et al.ll2009M ) 



IFGL J1826.2-1450 (LS 5039: lAbdo et al.ll2009of ) 
IFGL J2032.4+4057 fGvgX-3: IXbdo et al.ll2nn9j 



Formally, we associate five IFGL sources with low-mass X-ray binaries (LMXB). Three 
of them are also associated with globular clusters, and hence their combined emission from 
millisecond pulsars appears to be the more plausible counterpart of the IFGL sources (see 
§ I6.2.5P and we do not list LMXBs as a separate source class. One association corre- 
sponds to the Galactic Center source (cf. section I6.2.7p . And the remaining association, 
IFGL J1617.5-5105C (IE 161348-5055.1) is also associated with the SNR G332.4-00.4. 
Thus, none of the LMXB associations gives strong evidence that we indeed detect 7-ray 
emission from this source class. 



6.2.9. O stars, Wolf-Rayet stars, Luminous Blue Variable stars and open clusters 



The automated association procedure finds one O star (Cyg OB2-4) associated with 
IFGL J2032.2+4127, yet this source is know n to be a 7-ray pulsar (PSR J2032+4127, 



Abdo et all l2010ml . l2009cl : ICamilo et al. 



20091). In this c ase, the unusually large density 



of O stars in the Cyg 0B2 association jKnodlsedeil bood ) leads to source confusion, and 
ignoring this association we do not find any evidence for 7-ray emission from O stars in the 
IFGL catalog. 

Formally, two IFGL sources were associated with Wolf-Rayet stars. The first is the 
Galactic center source IFGL J1745.6— 2900c, which, as stated above, has many possible 
alternative associations, so we definitely cannot establish a physical link between the LAT 
source and t he Wolf-Rayet stars . The second is IFGL J2032. 4+4057 which has been identified 
as Cyg X-3 (lAbdo et al.ll2009q[ ). Cyg X-3 is a compact binary system that indeed comprises 
a Wolf-Rayet star, hence here the association is indee d correct. It is unli kely that the 7-ray 
emission is indeed arising from the Wolf-Rayet star (lAbdo et al. l2009q[ l. Thus we also do 
not find any evidence for 7-ray emission from Wolf-Rayet stars in the IFGL catalog. 

Two IFGL sources were associated with Luminous Blue Variable (LBV) stars: IFGL J1745.6— 2900c, 
the Galactic center source, which for the same reason as given above we do not consider as 
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a relevant association, and IFGL J1746.4— 2849c, which is associated with FMM 362 and 
the Pistol Star. However, the latter source is also associated with the PWN GO. 13— 0.11, so 
also here it is difficult to establish a physical link. 

Finally, we note that the LBV star rj Carinae was not formally associated by our proce- 
dure to the nearby source IFGL J1045.2— 5942. The formal posterior association probability 
for 7] Carinae is 0.76, hence below our threshold (0.8) for listing an association. The angular 
separation is 1.7' which is slightly larger than the 95% containment radius of 1.4'. Thus t] 
Carinae falls just outside the 95% error radius of IFGL J1045.2— 5942. On the other hand, 
IFGL J1045.2— 5942 was associated with the open cluster Trumpler 16 (the only association 
to an open cluster for the IFGL catalog), which besides rj Carinae houses many massive stars, 
similar to the Cyg 0B2 association. Recently, the young energetic pulsar PSR J2032+4127 
was found by the LAT in Cyg 0B2 (see above) and it is possible that young energetic pulsars 
are also hidden in Trumpler 16. 



6.2.10. Associations with Bright Source List sources 

Table [9] reports 195 associations between the 205 BSL (OFGL) sources and those of the 
IFGL catalog. Strictly applied, the automated association procedure found 186 associations. 
We have manually adopted nine more for BSL sources confidently associated with blazars 
or pulsars. Regarding the latter, the IFGL positions are formally 0, indicating that the 
source positions were assigned to the more-accurately known positions of the pulsars (see 
§111), affecting the statistical relationship between the OFGL and IFGL positions for these. 

The 10 BSL sources that do not have clear counterparts in the IFGL catalog are listed 
in Table [TOl Each is in the Galactic ridge, where for the IFGL catalog analysis we have 
recognized difficulties detecting and characterizing sources (§ 14. 7p . 



6.2.11. Associations with EGRET and AGILE sources 



The sources in the IFG L catalog have positional matches with 107 of the 271 3EG 
sources ( iHartman et al.lll999l ) (f our 3EG sources are each res olved into two IFGL sources) 
and 66 of the 188 EGR sources (ICasandjian fc Grenierll2008l ). A few more of the EGRET 
sources are close to, but not formally consistent with, IFGL source locations; the EGRET 
positions for the bright pulsars were offset from the true positions, for example. Almost 
all of the AGNs labeled in the 3EG catalog as good candidates are seen by LAT. One of 
the exceptions is 3EG J1230— 0247 (EGRc J1233— 0318), which was seen only early in the 
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CGRO mission. 

The IFGL catalog clearly does not account for a large fraction of the sources seen 
by EGRET. In light of the high sensitivity of the LAT and the fact that the LAT sees 
most of the EGRET catalog AGNs, which are known to be variable, the absence of more 
EGRET sources from the LAT catalog cannot be attributed primarily to time variability. A 
more likely explanation would seem to be that many of the EGRET sources were not discrete 
sources but were degree scale or larger diffuse structures not included in the model of Galactic 
diffuse emission used for analysis of the EGRET data. Its improved angular resolution and 
high photon statistics at GeV energies make the LAT far less sensitive to such structures, 
and the model of Galactic interstellar diffuse emission itself has incorporated far more detail 
than was available in the EGRET era. See § 14. 7^ however, for a discussion of how even some 
IFGL sources may be affected by the modeling of the diffuse 7-ray emission. 



All 47 of the sources in the first AGILE catalog (jPittori et al.ll2009l ) have corresponding 
sources in the IFGL catalog. A number of the lAGL sources map to multiple IFGL sources, 
and a few of the lAGL sources are close but not formally consistent in position with the 
IFGL sources. Nevertheless, the two present high-energy 7-ray telescopes do appear to be 
consistent in their detections of bright sources. 



6.3. Firm Identifications 



Firm identifications, indicated in the main table by capitals in the Class column, require 
more than a high-probability positional association. The strongest test for identification is 
time variability, either periodicity or correlation with variability seen at another wavelength. 
The 56 pulsars that have class PSR all show high- confidence (statistical probabili ty of chance 
occurre n ce less than 10~^) periodicity caused by the rota tion of the neutron star (lAbdo et al. 
2010ml B Is az Parkinson et al.ll2010l : lAbdo et al.ll2010bl ). Similar confidence levels apply to 
the t hree X-ray binary s ystems w hose orbital periods are detected in the LAT data: LSI -|-6 1 
303 kbdo et al.ll2009kl ). LS5039 JAbdo et al.ll2009oh . and Cygnus X-3 JAbdo et al.ll2009qh . 



With the large number of blazars detected by the LAT and the significant variability seen 
in many of these, the search for correlated variability that can provide firm identifications is 
a major effort that has not yet been carried out systematically for the LAT data. We have 
therefore chosen to list as firm identifications onl y those blazars for w hich publ ications exist 
showing such va riability. These are just 3C 273 (lAbdo et al.ll2010g]). 3 C 279 (llafrate et al. 
2009h. 3C 454.3 jAbdo et al.l l2009eh. PKS 1502+106 JAbdo et al.ll201ol ) and PKS 1510-08 



[e.g., 



Tramacerd l2008l ) . Additional studies will undoubtedly expand this hst. 
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Another approach to firm identification, slightly less robust than time variabihty, is 
morphology: finding spatial extent in a 7-ray source that matches resolved emission at other 
wavelengths. Some SNRs have measurable spati al extents in the L AT data and can be 



consid ered firm ide ntifications; here w e cite W44 (lAbdo et al.ll2010h[ ). W51C f lAbdo et al. 



2009jl ). and IC 443 (lAbdo et al.ll2010il ). The analyses for the IFGL catalog assume point-like 
emission, and so the positions and fluxes are not as well characterized as they would be in 
analyses that take into account the finite angular extents of these sources. A special case is 
IFGL J1322.0— 4515, which appears to b e part of one of the lobes of the emission from radio 
galaxy Centaurus A ( lAbdo et al.ll2010d[ ). Studies of source morphology are ongoing. 



The Large Magellanic Cloud (LMC) is an extreme example ( lAbdo et al.ll2009fi ) of the 
issue of analyzing an extended source with tools designed to find point sources. As described 
in § 16.2.41 the catalog contains 5 sources that are likely to be related to the diffuse 7-ray 
emission of the LMC. The LMC can be considered a firmly identified LAT source, but it is 
not a point source or an ensemble of point sources. 



TeV Source Associations 



IFGL sources that are positionally associated with sources seen by the ground-based 
TeV telescopes are of particular interest because the TeV band overlaps with the LAT energy 
range, suggesting the potential for such sources to be physically related. As described in Ta- 
ble [H we investigated associations with the sources in the TeVCat compilation of detections 
from ground-based observatories. The compilation is growing with time, and information 
about the sources is subject to updates and refinements, but at any given time TeVCat 
represents a snapshot of current knowledge of the TeV sky. 

The association analysis was made separately for extended (> 40' diameter) TeV sources 
(11) and smaller sources. The latter are typically much smaller than 40' and were treated 
like point sources for the association analysis. The 'TeV column of Table [2] lists associations 
with extended sources as 'E' and smaller angular-size sources as 'P'. As the table indicates, 
about half of the TeVCat sources have positionally-plausible associations with IFGL sources. 

The associatio ns include the Mila gro sources MGRO J1908+06, MGRO J20 19+37, and 
MGRO J2031+41 dAbdo et al.lboOTh, as wel l as MGRO J0632+17 and MGRO J2228+61 
( lAbdo et al.l l2009pl : iGoodman fc SinnisI 120091 ) , reported as part of a se arch for spatial cor - 
relations between the Milagro skymap and sources in the Fermi BSL (lAbdo et al.ll2009pl ). 
A number of the Milagro detections are pulsars with PWNs. The association between GeV 
7-ray PSRs and the PWNs visible in TeV 7-rays seems well established, as has already been 
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discussed in lAbdo et al.l (l2010inl ) . 



Other source classes are represented among the TeVCat associations. The LAT PSR 
IFGL J1023.0-5746 is spatially consistent with HESS J1023-575, itself not yet firmly iden- 
tified, but noted for its possi ble connection to the yo ung stellar cluster Westerlund 2 in the 
star-forming region RCW49 (lAharonian et al.l 120071 ) . Blazars, particu larly BL Lac object s 
are also solidly connected bet ween the GeV and TeV energy ranges (lAbdo et al.l l2009ml ). 
as are two starburs t galaxies (lAbdo et al.l l2010al). an d the HMXB sources LS I -|-61 303 



(lAbdo et al.ll2009kh and LS 5039 flAbdo et al.ll2009oh . SNRs such as W51C f lAbdo et al. 
2009jl ) and Gas A (lAbdo et al.l 12010 j ) also connect the two energy regimes. 



Of course sources that are positionally consistent between the LAT and TeV telescopes 
but have no obvious associations with objects at longer wavelengths are also of interest; 
Table |2] lists a few of these. Establishing a physical connection through spectral or variability 
studies may help determine the nature of these sources. 



As discussed in Section I6.2.7[ the Galactic Genter region is particularly complex. In- 
vestigations of_Jh^_j^socw^ions with the TeV 7-r ay sources known, in th is region - HESS 
J1745-290 f|Acero et al.l l2009bh. HESS J1745-303 jAharonian et al.lboosi ) and HESS J1741-302 



( iTiboUa et al.ll2008l ) are outside the scope of this paper and will be discussed elsewhere. 



8. Conclusion 

The 1451 sources in this First Fermi-LAT catalog (IFGL) represent the most complete 
understanding to date of sources in the GeV sky. The catalog clearly contains a number of 
populations of 7-ray emitters. It offers a multitude of opportunities for additional research, 
both on sources with likely associations or identifications and on those sources that remain 
without apparent counterparts. 
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A. Estimation of Detection Threshold 

An approximate, but reasonably accurate, expression for the detection threshold for a 
point source at any point in the sky can be obtained by assuming that the diffuse back- 
ground is locally uniform and considering only one source. We start by constructing the log 
Likelihood function given by a sum of the logarithms of Poisson probabilities for detecting 
Hi events in some bin i. where i labels position in the sky and energy, when the model 
predicts A^: logL = '^■riilog Xi — Aj. In the above approximation we can replace the sum 
over i by an integral over energy and the angular separation between the source location 
and the event direction, {E,9). The expectation value for the detected counts density is 
n{E, 0) = ToAeff(E) [S{E)PSF{e, E) + B{E)] where S{E) and B{E) describe the spectra of 
the source and background events, PSF(6', E) the spatial distribution of events from a point 
source, Aes{E) the effective area of the instrument and Tq the observation duration. To 
estimate TS for the detection in this case we calculate the Likelihood under two hypotheses: 
the maximum likelihood hypothesis, L*, where the source term is included in the model, and 
the null hypothesis, where the counts are presumed to arise from the background only 
(also assuming the number of source counts is small over the full ROI in comparison to the 
total background counts). From these we derive TS = 2(logL* — logL°). Writing the local 
source-to-background ratio as g{9,E) = S{E)PSF{9, E)/B{E) we get: 

/■-Emax / /"I" \ 

TS = 2To / A,s{E)dEl B{E) [1 + g{e, E)]log[l + g{e, E)] dQ - S{E){Kl) 



This preprint was prepared with the A AS macros v5.2. 
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/•log Emax 

= To / W{E)d\ogE (A2) 

= 2EA,s{E)BiE) [ [l + g{e,E)]hg[l + g{e,E)]~g{e,E)dn (A3) 

Here is the contribution to TS per unit logi?. It is illustrated in Figure [T8] 

(dashed line) for a power-law source spectrum with index 2.2. At low energy (below 1 GeV) 
faint sources are always background limited, i.e., g{d, E) is small even at = 0. In that limit, 
the dominant term of the integrand in the integral over 6 in Eq. IA3I is g{6, Ey /2. Noting 
that PSF{e, E fdVt = Csh(^)cx(^)~^ in which a{E) is the angular resolution and Csh(^) 
is a shape factor weakly dependent on energy, we get 

This explicitly shows that the weight is proportional to the ratio of source counts over 
background counts within the angular resolution. The strong improvement of the PSF with 
energy (approximately as E~^'^) means that for a E'"^'"^ source (which is not far from the 
background spectral shape) W{E) oc Ae{i{E)E^-^. Since A^s^E) improves with energy up to 
1 GeV, this explains the rising part of W{E). At high energy (above several GeV) the PSF is 
narrow enough that even the faint sources are limited by their own count rates (^'(0, E) > 1) 
and Eq. IA4l no longer applies. When the source density is large, a first-order way to account 
for confusion is to limit the integral over angles in Eq. lASI to Omux such that tcO'^^^ = ^^tot/^src 
is the average solid angle per source. This is shown by the solid line on Figure [181 The effect 
is of course larger where the PSF is broader, at low energy. 



Setting TS = 25 in Eq. lAlt assuming a given source spectral shape and solving for 
the normalization of the source spectrum provides the detection threshold. The spatial 
dependence of that threshold is shown in Figure [191 The Galactic diffuse and isotropic 
backgrounds are taken from the model (§ [3]). Tq for each point is derived from the Fermi 
pointing history during the 11 months, and the source spectrum is assumed to be E'"^'"^, the 
average spectral slope of the IFGL sources. Although the nonuniform exposure affects this 
map somewhat, the dominant factor is the strong diffuse emission along the Galactic plane. 
Because of the strong energy dependence of the PSF, the detection threshold depends very 
sensitively on the spectral index as well. Figure [201 illustrates this in terms of the photon 
flux above 100 MeV, which ranges from 10^^ to 4 x 10^*^ ph cm^^ s^^ going from very hard 
(r = 1.5) to soft (r = 3) sources. The recipe for source confusion, given above, is used. We 
note that a few sources are below the line. This can happen for purely statistical reasons, 
or because the background and exposure depend a little on the direction, even after taking 
out the Galactic plane. 
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Table 10. Bright Source List (OFGL) Sources Without IFGL Counterparts 



Name OFGL 


I 


b 


J1106.4-6055 


290.52° 


-0.60° 


Ti 1 1 ^ s fil ns 


9Q1 


u.oo 


J1615.6-5049 


332.35 


-0.01 


J1622.4-4945 


333.87 


-0.01 


J1634.9-4737 


336.84 


-0.03 


J1741.4-3046 


357.96 


-0.19 


J1821.4-1444 


16.43 


-0.22 


J1836. 1-0727 


24.56 


-0.03 


J1900.0+0356 


37.42 


-0.11 


J1024.0-5754 


284.35 


-0.45 



0.8 
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Fig. 18. — Theoretical contribution {W{E) of Eq. IA3I1 to Test Statistic per Ms and per 
\og{E) interval as a function of energy for a power-law source over the average background 
at |6| > 10°. The assumed photon spectral index is 2.2. The dashed line is for an isolated 
source. The full line includes approximately the effect of source confusion. 
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Fig. 19.— Flux (E>100 MeV in ph cm-^ s'^) needed to reach TS = 25 in the LAT data 
for the 11-month time range considered in this paper, as a function of position in Galactic 
coordinates. The assumed photon spectral index is 2.2. 




Spectral index 



Fig. 20. — Photon flux above 100 MeV of sources at |6| > 10° as a function of spectral index. 
The dashed line shows the theoretical detection threshold at TS = 25 for an isolated source 
over the average background at |6| > 10°. The full line includes approximately the effect of 
source confusion as on Figure [181 
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B. Association Method 

The method implemente d for automat i c asso ciation of the IFGL sources essentially 



follows the ideas developed by iMattox et al.l (Il997l ) for the identification of EGRET sources 



with fiat-spectrum radio sources. It makes use of Bayes' theorem to compute the posterior 
probability Pjfc(A|r, 0) that a counterpart i from a list of potential counterparts supplied in 
the form of a counterpart catalog is the correct association of a LAT source k: 

' '"^^ p,fc(r,0|A)P,(A)+p,fc(r,0|A)Pi(A)' ^ ^ 

Pi{A) is the prior probability that counterpart i is detectable by the LAT, Pi(A) = 1 — 
Pj(A) is the prior probability that counterpart i is not detectable by the LAT, pik{r,(j)\A) 
is the probability density for the detectable counterpart i to have an angular separation r 
and position angle from a LAT source k, and pik{r,(j)\A) is the probability density for 
source k being only by chance situated at an angular separation r and position angle (p from 
counterpart i. 

Under the assumption that source location confidence region is centrally peaked, the 
probability density pik{r,(j)\A) is obtained by differentiation of the probability Pife(r, 0|A) 
that the LAT source k is located at an angular separation smaller than r under the position 
angle 0: 

Pifc(r,0|A) = — — . (B2) 

rdrdcp 

By defining 

Afc = Inmax(Lfe) - lnLfc(r, 0) (B3) 

as the difference between the log-likelihood maximum of the LAT source k and the log- 
likelihood at position (r, 0) and by maki ng use of Wilks's theorem that 2Afc is distributed 



as xi '^ull hypothesis IWilkd (Il938[ ) one can write 



oo 

P,fc(r,0|A) = l- / xl{x)dx = l-e-''K (B4) 



2At 



We approximate as an elliptical paraboloid which is defined by the 95% confidence 
elliptical error region: 



/cos2(0 - 0fc) sm^{(p-(pk\ 



A. = - ( + ] (B5) 



2 V al hi 



where and hk are the axes of the ellipse at 1 a, smaller than the semimajor (Conf _95_SemiMaj or 
in the FITS version of the IFGL catalog; App. D) and semiminor (Conf _95_SeiniMinor) axes 
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at 95% confidence by a factor a/— 2 log(0.05) = 2.45, is the position angle of the error 
elhpse (Conf _95_PosAng) and (p is the given position angle between LAT source k and the 
counterpart i. This can be transformed to 

The chance coincidence probability density Pifc(r, 0|A) is determined from the local den - 



sity Pk of counterparts around LAT source k as proposed by [Sutherland fc Saunders! (119921 ) : 



To compute this density we count the number of counterparts A^^ in the counterpart catalog 
under consideration within a radius of ro = 4° around the location of the LAT source k and 
divide by the solid angle of the search region: 

Pk = ^ (B8) 

Note that the counterpart i is included in A^^ which guarantees that A^jt > 1. 

As a last step, we implement the reasonable condition that a counterpart i cannot be 
associated with more than one LAT source. This is done by introducing A^lat + 1 mutually 
exclusive hypotheses (A^lat = 1451 being the number of sources in the first year catalog): 

Hfc: Object i is a counterpart of LAT source k and of none of the other LAT sources. 
H_: Object i is not a counterpart of any LAT source. 

The probabilities for these new hypotheses are computed using 

Pifc(Hfe|r,0) = P,fc(A|r,0) J]P,fc,(A|r,0) 

k'^k 

P,(H_|r,0) = J]P,fc,(A|r,0). (B9) 

k' 

where 

P,fc(A|r,0) = l-P,fe(A|r,(/)). (BIO) 

Since we dropped from the set of hypotheses all of the cases where an object i is associated 
to more than a single LAT source, the sum over all probabilities 



S^ = 5^P.fc(Hfc|r,0) + P,(H_|r,0) (Ell) 
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is < 1, and we thus renormalize using 

Pifc(Hfc|r,0) = ^ (B12) 

to obtain the posterior probabihty that object i is a counterpart of the LAT source k and of 
none of the other LAT sources. Practically, Si < 1 only if the error ellipses of neighboring 
LAT sources overlap, which is rather unlikely. Thus, to a good approximation we have 
Pjfc(Hfc|r, 0) = Pjfc(A|r, (^) for basically all sources. 

The above procedure that leads to the computation of Pjfc(Hfc|r, 0) for a specific catalog 
of counterpart candidates has been implemented in the ScienceTools executable gtsrcid. 
We used version v2r2p3 of this tool for counterpart association for the first year catalog. To 
simplify, in what follows we will write Pik instead of Pik(Rk\r, 0) for the posterior association 
probability of LAT source k with object i of a given counterpart catalog. 



C. Calibration of prior probabilities 

Before Eq. (IBip can be used for source association, prior probabilities Pi{A) have to 
be specified for each counterpart i. Here we make the simplifying assumption that within a 
given counterpart catalog the prior probabilities for all sources i are identical: 

P,(A) = P(A). (CI) 

To assign P(A) for a given counterpart catalog, we require the relation 

A^false= (C2) 

to hold, where Nfs,\se is the number of false associations that have posterior probabilities 
Pik above our selected threshold Pthr = 0.8. For this purpose we determined by means of 
Monte Carlo simulations of 100 fake LAT catalogs the expected number of false associations, 
(iVfaise), as function of P(A). For a given P(A), we obtained (iVfaise) by counting the number 
of sources in each fake catalog that have been associated, and by dividing this number by 
100, i.e. by the number of fake catalogs that have been simulated. P(A) has then been 
varied until A^faise = (-^faise) was fulfilled, which then fixed the proper prior probability P(A) 
for the counterpart catalog. 

The fake catalogs were created by randomly displacing IFGL sources within a ring from 
2° to 10° in radius around their nominal position. Since the IFGL catalog comprises a 
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distinct population of Galactic sources that obey a rather narrow latitude distribution, we 
limited source displacement in Galactic latitude to b ± 6max; where 

&max = '^max (^1-0 " SCCh^ (^y^ ^ , (C3) 

'"max = 10°, b is the Galactic latitude of the IFGL source, and 60 = 5° is the angular scale 
height above the Galactic plane for which the latitude displacement is reduced. We further 
required 6max > 0.2° to allow for a non-zero latitude displacement of sources in the Galactic 
plane, and required any source to be shifted by at least r^in = 2° away from its original 
location. For illustration, we show in Figure [21] the locations of the sources in the 100 fake 
catalogs that were used for calibration. 

As example, we show A'faise (determined from Eq. IG2p and (iVfaise) (determined from the 
simulation) as function of -P(A) for the CRATES catalog of flat spectrum radio sources in 
the left panel of Figure [221 The intersection of both curves determines the prior probability, 
which in this case has been determined to -P(A) = 0.33. We also show in Figure [22] the 
number of associations that is computed using 

as solid black lines. 

Although the calibration has been performed for Pthr = 0.8, it turns out that, once 
the prior probabilities are set, Eq. (]C2p is fulfilled for a large range of posterior probability 
thresholds. We illustrate this property in the right panel of Figure [221 which shows for the 
CRATES catalog (iVfaise) and A/'faise as function of Pthr for P(A) = 0.33. Obviously, both 
quantities are in good agreement. 



D. Description of the FITS version of the IFGL catalog 

The FITS format version of the IFGL cataloj^ has three binary table extensions. The ex- 
tension LAT_Point_Source_Catalog Extension has all of the information about the sources, 
including the monthly light curves (Tab. [TTl) . 

The extension Hist_Start lists the Mission Elapsed Time (seconds since 00:00 UTC on 
2000 January 1) of the start of each bin of the monthly light curves. The final entry is the 
ending time of the last bin. 



The file is available from the Fermi Science Support Center, 



http://fermi.gsfc.nasa.gov/ssc 




Fig. 21. — Locations of sources in the 100 fake IFGL catalogs used for calibration of prior 
probabilities. 





Posterior Probability Threshold 



Fig. 22. — Number of false associations as function of prior probability -P(A) for Pthv = 0.8 
(left panel), and as function of posterior probability for P(A) = 0.33. Dashed red lines 
correspond to (iVfaise), solid blue lines represent A^faise (left axis). The solid black lines show 
A^ass (right axis). 
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The extension GTI is a standard Good-Time Interval listing the precise time intervals 
(start and stop in MET) included in the data analysis. The number of intervals is fairly 
large because on most orbits (~95 min) Fermi passes through the South Atlantic Anomaly 
(SAA), and science data taking is stopped during these times. In addition, data taking is 
briefly interrupted on each non-SAA-crossing orbit, as Fermi crosses the ascending node. 
Filtering of time intervals with large rocking angles, other data gaps, or operation in non- 
standard configurations introduces some more entries. The GTI is provided for reference 
and would be useful, e.g., for reconstructing the precise data set that was used for the IFGL 
analysis. 
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Table 11. LAT IFGL FITS format: LAT_Point_Source_Catalog Extension 



Column 


Format 






Unit 


Description 




18 A 












RA 


E 


deg 








Riglit Ascension 


DEC 


E 


deg 








Declination 


GLON 


E 


deg 








Galactic Longitude 


GLAT 


E 


deg 








Cfilactic TjfititiiHp 


Grtnf fiS SpmiA/Ta inr 


E 


deg 








TiOTic radin*? nx prT'or p11iti*5p at fiSvri ponfidpncp 


Conf_68_SemiMinor 


E 


deg 








Sfiort rfl/1in<s OT prroT* pllin*ip at flSvfi cnnfidpnpp 


Conf_68_PosAng 


E 


deg 








Position angle of the 68% long axis from celestial North, 

positive toward increasing RA (eastward) 


Conf 95 SciiiiMcijor 


E 


deg 








Long radius of error ellipse at 95% confidence 


Conf_95_SciiiiMiiior 


E 


deg 








Short radius of error ellipse at 95% confidence 


Conf Q5 PosAne- 


E 


deg 








PriGi + ir^n civicrlp nT trip Q^^Vi Inner dYiG TVi~ivn ^plpG^'icil IXTrii'+n 
X ^^DlUHJll dllglC Ol liiiC i70 /U lOii^ lHJlll ^jClCOUldl i^H^lUilj 

T\^^*;^^■^T^p ^'mjimyn inf^ypadno" R A ( pa*;^■^]^^^^T■H i 




E 










Source significance in sigma units (derived from Test Statistic) 


Pivot Energy 


E 


McV 






Energy at which error on differential flux is minimal 


Flux Density 


E 


cm 


-2 


MeV-1 s-l 


Differential flux at Pivot Energy 


Unc_Flux_DGnsity 


E 


cm~ 


-2 


MeV-i s-i 


1 (T error on differential flux at Pivot.Energy 


SriPfi^Tfi 1 Tn H PY 


E 










Rpct fit nrnjirpi" Imjtr dnr^p 


T Tup SnpptTal TrirJpv 


E 










1 (T prror on lip*it fit tiowpt" law siloTip 


FluxlOOO 


E 


cm~ 


-2 


s" 


-1 


Integral flux from 1 to 100 GeV 


Unc_FluxlOOO 


E 


cm" 


-2 




-1 


1 (7 error on iutegrctl flux from 1 to 100 GeV^ 


Ensrgy Flux 


E 


erg cm 




Energy flux from 100 MeV to 100 GeV 


Unc Energy Flux 


E 


erg cm" 




1 (7 error on energy flux from 100 ]VIeV to 100 GeV 


(111 wfltiivp TnHpy 


E 










IVTpa^iiTP of liow ^npftniTTi follnws nrtwpv-law ^'fnrrpntlv cjirnrilp 'v^ 


FlioxSO-lOO 


E 


cm 


-2 


s" 


-1 


Integral flux from 30 to 100 MeV (not filled) 


Unc_Flux30_100 


E 


cm~ 


-2 


s" 


-1 


1 o" error on integral flux from 30 to 100 Iv^eV^ (not filled) 


Sqrt_TS30_100 


E 










Square root of the Test Statistic between 30 and 100 MeV (not filled) 


Fluxl00_300 


E 


cm" 


-2 


s" 


-1 


Integral flux from 100 to 300 MeV 


TInr Fliixl nfl SflO 

KJ 11L^_-L 1 Ll^V ±^ W VJ • f ^ J \J 


E 


cm" 


-2 


s" 


-1 


_L (_/ KzL L^^L Wli lllLrs_.tlI- <X1 11 Ll^V llWlll -L \J\J LLJ rj\J\J IVXO V 


Sqrt_TS100_300 


E 










Square root of the Test Statistic between 100 and 300 MeV 


Flux300.1000 


E 


cm~ 


-2 


s" 


-1 


Integral flux from 300 MeV to 1 GeV 


Unc.Flux300.1000 


E 


cm~ 


-2 


s' 


-1 


1 <T error on integral flux from 300 MeV to 1 GeV 


Sqrt.TS300.1000 


E 










Square root of the Test Statistic between 300 MeV and 1 GeV 


FluxlOOO.3000 


E 


cm" 


-2 


s" 


-1 


Integral flux from 1 to 3 GeV 


Unc.Fluxl000.3000 


E 


cm~ 


-2 


s' 


-1 


1 cr error on integral flux from 1 to 3 GeV''' 


Sqrt.TS1000.3000 


E 










Square root of the Test Statistic between 1 and 3 GeV 


FluxSOOO.lOOOO 


E 


cm" 


-2 


s" 


-1 


Integral flux from 3 to 10 GeV 


Unc.Flux3000.10000 


E 


cm~ 


-2 


s" 


-1 


1 cr error on integral flux from 3 to 10 GeV 


Sqrt.TS3000.10000 


E 










Square root of the Test Statistic between 3 and 10 GeV 


FluxlOOOO.lOOOOO 


E 


cm~ 


-2 


s" 


-1 


Integral flux from 10 to 100 GeV 


UncFluxlOOOO.lOOOOO 


E 


cm~ 


-2 


s" 


-1 


1 a error on integral flux from 10 to 100 GeV^ 


Sqrt.TSlOOOO.lOOOOO 


E 










Square root of the Test Statistic between 10 and 100 GeV 


Variability Jndex 


E 










Measure of source variability (currently simple x^) 


SigniLPeak 


E 










Source significance in peak interval in cr units 


Flux.Peak 


E 


cm~ 


-2 


s" 


-1 


Peak integral flux from 100 MeV to 100 GeV 


Unc.Flux.PeaJj 


E 


cm~ 


-2 


s" 


-1 


1 cr error on peaJj integral flux 
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Table 11 — Continued 



Column 


Format 


Unit 


Description 


Time_Peak 


D 


s (MET) 


Time of center of interval in which peak flux was measured 


Peakjnterval 


E 


s 


Length of interval in which peak flux was measured 


Flux_History 


HE 


— 2 —1 

cm s 


Integral flux from 100 MeV to 100 GeV in each interval 


Unc_Flux_History 


HE 


— 2 —1 

cm s 


Error on integral flux in each interval using method 

indicated in Unc_Flag_History column and added in quadrature 

with 3% systematic component. 


Unc_Flag_History 


IIB 




1 if it is half of the difference between the 2a upper limit 

and the maximum-likelihood value given in Flux_History, if it is the 

Ict uncertainty derived from a signiflcant detection in the interval 


OFGL_Name 


18A 




Name of corresponding OFGL source, if any 


ASSOC.GAMl 


ISA 




Name of likely corresponding lAGL source 


ASSOC.GAM2 


ISA 




Name of likely corresponding 3EG source 


ASSOC.GAM3 


ISA 




Name of likely corresponding EGR source 


TEVCAT.FLAG 


A 




P if positional association with <40' source in TeVCat 

E if associated with a more extended source in TeVCat, N if no TeV association 


CLASSl 


3A 




Class designation for associated source; see Table |6] 


CLASS2 


3A 




Second class designation for associated source 


ASSOCl 


24A 




Name of identified or likely associated source 


ASSOC2 


24A 




Alternate name of identified or likely associated source 


Flags 


I 




Source flags (binary coding as in Table |4]l 



''The upper limit is set equal to if the flux in the corresponding energy band is an upper limit (TS < 10 in that band). 
The upper limits are 2 a. 



